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ABSTRACT 


The advent of the space shuttle (STS) will make it possible 
to carry out relativistic experiments in a laboratory in free 
fall. One of the most important such experiments , because it 
is essential to the foundations of general relativity, is the 
test of the universality of the ratio of passive gravitational 
to inertial mass for bodies, which is known as the EOtvOs ex- 
periment. This report presents analyses of a direct force- 
balance technique for carrying out this experiment in space, 
which is intended to give sufficient sensitivity to allow in- 
vestigation of the gravitational interactions of energy stored 
in the weak interaction. It is found that a sufficiently sensi- 
tive experiment may be possible in which the apparatus is al- 
lowed to float in the payload bay of the shuttle, although maxi- 
mum sensitivity requires a fully autonomous, free-flying experi- 
ment. 


The heart of this experiment is an exceedingly sensitive dual 
accelerometer containing two proof masses constructed of the 
materials whose EOtvOs ratio is to be compared. It is hoped to 
use an electrostatically-suspended accelerometer for this appli- 
cation; but it is impossible to test such an accelerometer in a 
terrestrial laboratory if electrostatic forces are used to sup- 
port the proof mass, because of cross-coupling between the support 
and sensitive axes. For use in development of the accelerometer, 
a magnetic microbalance is proposed in which the weight of the 
proof mass is supported by magnetic forces which vary very slowly 
with the proof mass position. It may then be possible to build 
an electrostatic suspension to control the proof mass with elec- 
trode spacings, voltages, etc., similar to those which are needed 
in free fall. It is shown that at least two different mechaniza- 
tions of the magnetic suspension may be feasible. 
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CHAPTER ONE 

THE EOTOVOS EXPERIMENT 

I. 1 Introduction 

In Newtonian gravitational theory, the force acting on a 
particle of inertial mass m^, due to the gravitational field of 
another particle of mass. ^ is given by 

F = m^r = - GM^M2 r~^r [1] 

where G is the gravitational constant, r is the 
position vector of the particle relative to the source of the 
field and M^ is the gravitational mass of the particle. As has 
been noted by Bondi^, three logically distinct concepts of mass 
occur in this equation: (i) inertial mass m^, which determines 

the acceleration of the particle under a given force; (ii) passive 
gravitational mass M^, which determines the force on the particle 
in a given gravitational field; and (iii) active gravitational 
mass M 2 , which determines the strength of the gravitational field 
generated by a particle. The equivalence of active and passive 
gravitational mass is a consequence of the law of action-reaction 
[which is not a Lorentz-invariant concept, so that it must be 
handled with care in relativistic generalizations of the theory] , 


but the universality of gravitational phenomena depends on the 
separate assumption that the Eotvos ratio 

k = M/m [2] 

is a constant for all bodies (at least, for all bodies 
which are so small that gravity-gradient effects in the external 
gravitational field are unimportant) . The value of G is normally 
so chosen that this ratio is unity. 

The hypothesis that k is independent of the chemical composition 
or internal structure of a sufficiently small body is a statement of 
the Weak Principle of Equivalence, which is an axiom in general 
relcitivity and in most other theories of gravitation, although in 
some (for example, the Brans-Dicke cosmology^) the ratio is allowed 
to vary from point to point in spacetime . It should be noted that 
the Weak Principle is necessary but not Sufficient for the validity 
of general relativity (and especially for its geometric interpretation 
in terms of the metric of spacetime) ; the theory developed by 
Einstein requires the Strong Principle of Equivalence, which states 
that no self-contained physical measurement, carried out in an 
infinitesimal region of spacetime, can distinguish between gravitation 
and a suitably chosen inertial acceleration.' As a corollary, a 
free-fall coordinate system, if sufficiently limited in extent, 
is physically absolutely equivalent to an inertial frame (local 
Lorentz '.frame)' . 
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The Strong Principle is subject to direct experimental proof 
only if it is believed that it is presently possible to enumerate 
every feasible type of physical measurement. However, it is 
possible that the Weak Principle implies the Strong Principle, 
a suggestion known as the Schiff Conjecture^ which has, as yet, 
been satisfactorily demonstrated only for restricted cases. ^ If 
it can be proven, the Weak Principle will become the primary 
foundation of gravitational theory. 

The Weak Principle of Equivalence is obviously supported by 
the common observation that all bodies fall with the same accelera- 
tion under gravitation (as in the famous and perhaps apocryphal 
test by Galileo at the Leaning Tower in Pisa) . This phenomenon 
was demonstrated dramatically on television by David Scott on 
Apollo 15, when he dropped a feather and a hammer together to 
the lunar surface. Accurate evidence for the validity of the 
principle however consists of experiments of the Ebtvos type, 
whose distinguishing characteristic is that they are null measure- 
ments, taking advantage of some situation in which there is a 
nominal balance between inertial and gravitational forces in order 
to achieve quite remarkable precision. Given two bodies, of 
materials A and B, these experiments may be regarded as measurements 
of the Eotvos number 

ri (A,B) = k (A) - K(B) [2A] 


Eotvos himself used the fact that, in a horizontal plane 
determined by a plumb bob, there is (in the northern hemisphere) 
a small southward component of centrifugal acceleration, due to 
the diurnal rotation of the Earth, and an equal but opposite 
component of gravitation. The magnitude of these horizontal 
accelerations reaches a maximum of about 1.7 cm/sec^ at 45° latitude 
A torsion pendulum was constructed in which bobs of different 
materials were attached to opposite ends of an arm which was 
suspended by a fiber at its balance point. If r) for the two bobs 
differed from zero, a small torque would be produced about the 
vertical, which would affect the rest orientation of the arm 
relative to the laboratory. Because this Eotvos torque was a 
sinusoidal function of the azimuth angle of the arm, the arm 
orientation would not then change by exactly 180° when the fiber 
suspension was rotated through precisely that angle. In a series 
of painstaking experiments between 1889 and 1922, Eotvos showed^ 
that, for a variety of pairs of materials, g differs from zero by 
at most a few parts in 10^, 

. C , 

In 1962, Dicke et al D built an improved torsion balance, which 
possessed triangular symmetry to reduce the effects of local gravity 
gradients and which was designed to take advantage of the balance 
which exists between the gravitational pull of the sun and the 
centrifugal force due to the motion of the Earth in its orbit. 

The accelerations involved here were somewhat smaller (c. 0.6cm/sec2 
than in the experiments of Eotvos, but they were modulated with a 


24-hour period by the rotation of the Earth, allowing frequency 
discrimination against some interfering torques and, more importantly 
eliminating the need to disturb the system by rotating the suspension 
relative to the laboratory. The accuracy achieved was about one 
part in 10^, in a comparison of gold and aluminum test masses. 

Using a modified form of the Dicke experiment, Braginski? in 
1972 improved the accuracy further-, by about one order of magnitude. 

As discussed in Chapter II, the advent of space technology 
has opened up opportunities for Eotvos experiments of unprecedented 
accuracy. 

1.2 Possible Violations of the Equivalence Principle 

General relativity is a theory of gravitation which is 
self-consistent, complete, and in agreement with all experiments 
to date. Moreover, it is compellingly elegant, in a way unmatched 
by any other theory which has been proposed. Even if other gravita- 
tional theories are to be considered, there are strong arguments^ 
that the only viable theories are metric theories — i.e., theories 
in which spacetime possesses a metric which satisfies the Equivalence 
Principle by exhibiting locally Lorentz frames. The best confirma- 
tion of the Equivalence Principle (and especially of the Strong 
Principle) lies in the nature of the theories to which it leads, 
rather than in the direct experimental evidence for it. Nevertheless 
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an unsatisfactory aspect of metric theories of gravitation is that, 
despite much effort which has gone into unified field theories, 
gravitation remains conceptually quite distinct from the other fields 
of physics, being regarded as an expression of the curvature of 
spacetime. Furthermore, as will be shown, there are difficulties 
with the Equivalence Principle itself which suggest that anomalies 
might be revealed by continued experimentation. 

In designing such tests, a theoretical framework would obviously 
be useful which predicted where violations might be found. It is 
not reasonable to discuss tests of the Equivalence Principle in 
terms of a theory which assumes its validity. What is needed is 
a gravitational theory which starts with minimal assumptions 
(e.g., Lorentz covariance and reduction to Newtonian theory in 
that limit) , identifies clearly the additional assumptions which 
lead to a metric theory, and then systematically explores the 
consequences of changing those assumptions. The spin-two field 
theory of gravitation 9 ' 10 , which leads to general relativity when 
non-linear terms due to the energy density of the gravitational 
field itself are included, may provide a model for such a development. 11 

In the absence of a comprehensive theory, the best that can 
be done is to consider some isolated possible anomalies. 

1.2.1 Differences in Composition: The Weak Interaction 

The classical Eotvos experiment is a test of the possibility 
that the ratio of passive to inert mass of a body depends on its ■ ' 
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STRONG INTERACTION 




ELECTROMAGNET I C 
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chemical composition. In more fundamental terms, Jie hypothesis 
to be tested is that one or more of the forms of energy which make 
up matter exhibits an anomalous value of this ratio. Examination 
of this hypothesis will show the motivation for further improvements 
in the accuracy of the experiment, as well as providing some important 
design considerations. 

Ordinary matter may be regarded as consisting of energy stored 
in a combination of the following forms: 

(i) The elementary particles (both real and virtual) which 

make up atoms, particularly protons, neutrons and electrons. 

(ii) The strong interaction, which binds protons and neutrons 
together to form nucleii. 

(iii) The electromagnetic interaction, which binds electrons 
and nucleii together to make atoms, atoms together to 
make molecules and crystal lattices, and which weakens 
the strong attraction of protons in a nucleus. 

(iv) The weak interaction, which is responsible for 8 -decay 
processes. 

(v) The gravitational interaction, which holds stars together 
and binds planets and stars into solar systems and solar 
systems into galaxies, but which is very weak indeed 
on the laboratory scale. 

Eotvos experiments using bodies of widely different composition 
allow conclusions to be drawn concerning the passive mass of some 
of these forms of energy (excluding the gravitational interaction,, 
which is discussed in Section 1.2.2). Energy conservation requires 
conservation of passive mass in reversible transformations between 
forms of energy, since otherwise it would in principle be possible 
to convert a system to its heavy form, extract energy by lowering 
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it in a gravitational field , and then lift it after conversion 
back to its light form, allowing a perpetual motion machine. However, 
this does not prohibit equal and opposite changes in passive mass 
amongst components of such transformations , so that anomalous 
results of the Eotvos experiment are possible. 

In order to estimate the accuracy required in such experiments, 
assume for simplicity that only the i th form of energy exhibits an 
Eotvos anomaly, with k^ the ratio of passive to ; inert mass. If 
(A) is the fractional contribution of this energy to the overall 
inert mass-energy of body A, Eq. (2A) reads 

n (A,B) = (1 - k ± ) [a ± (B) - at (A) ] [3] 

Since the mass fractions of the different energy forms 
add to unity for a given body, the difference in the particle mass 
fractions for two bodies is equal, apart from sign, to the difference 
in the binding energy mass fractions. Figure I shows the difference 
between the total binding energy for an element of atomic number Z 
and that of beryllium, which is arbitrarily chosen as a reference. 
This curve, which is obtained from the semi-empirical mass formula 
of Weizsacker 1 ^ , may be useful in designing tests of the hypothesis 
that particles, of whatever type, exhibit an anomalous Eotvbs ratio. 
FOr example, it is clear that this anomaly, if it exists, would 
b£ maximized in a test which compared beryllium with copper or 
iron. Such a choice would in fact give a test 6 times as sensitive 
as one using gold and aluminum, with no other changes in the 
experimental conditions . 
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The total binding energy is made up basically of the electromagnetic 
and strong interaction energies, the latter being negative since 
it is an attractive, force between nucleons. The weak and gravitational 
interactions are negligible by comparison. The electromagnetic 
and strong components of the difference in binding energy of an 
element of atomic number Z and that of beryllium are also shown 
in Fig. I, where it is seen that the optimum choice for test of 
the hypothesis that one of these energies has an anomalous Eotvos 
ratio is beryllium and an element of high atomic number, such as 
gold or uranium. The achievable improvement over gold and aluminum 
is by a factor of 2.4 for the strong interaction, but only by 
about 50% for the electromagnetic interaction. 

The shape of the curves in Fig. I suggest that an Eotvos 
experiment in which beryllium is used as a standard for comparison 
with both copper and gold would give optimum sensitivity fc_ each 
of the. three above forms of energy — and if an anomaly were 
detected, it would be possible to determine, from the relative 
magnitude of the results for copper and gold, which energy was 
responsible. 

The Dicke test, showing that n(Au,Al) < 10 , allows the 

following conclusions to be drawn from Fig. 'I: 

^particles i- s l ess than 2 x 10 "^ 
kgtrong l ess than 5 x 10 " 9 

^electromagnetic ^- s l ess than 6 x 10 ^ s 
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Note that this simple analysis applies to particles only when 
they are considered as a group. Drawing conclusions from the 
Eotvos experiment about the passive mass of particular types of 
particles (e.g. , neutrons) requires a more detailed analysis, as 
the differences in mass fractions are not necessarily given with 
sufficient accuracy by the semi-empirical mass formula and, in 
any case, depend on the particular isotopic composition of the 
samples tested. An analysis by Schiff-^ suggests that Eotvos 
experiments to date have demonstrated with convincing precision 
that protons, neutrons and electrons (and their antiparticles) 
do not exhibit anomalous gravitational behavior. 

Energy stored in the nucleus by virtue of the weak interaction 

■ 7 

is typically of order 10 times less than that due to the strong 
interaction-^; in other words, weak interaction energy may contribute 

9 

a fraction of order one part in 10 to the mass of an atom. In 
view of the uncertainties in this calculation, and of the fact that 
only differences in the mass fraction between different materials 
contribute in an Eotvos experiment, it is clear that such experiments 
to date are insufficient to allow any statement about the passive 
mass of weak energy. On the other hand, the Eotvos experiment in 
orbit which is discussed in Chapter II, which is intended to have 
a sensitivity of about one part in 10^, is fully capable of detecting 
any substantial violation of the Equivalence Principle by weak 
energy. It is this possibility which provides the strongest motivation 
for the Eotvos experiment in space. 
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1.2.2 The Gravitational Self-Energy 

For an approximately spherical body of radius r, the fractional 
contribution of the. energy in its gravitational field to the overall 
inert mass is of order 

a g = — 

rc 2 [4] 

where c is the speed of light. This number is of order 

— 25 

10 for laboratory bodies, much too small to be of any experimental 
interest, but it is of order 10“ 9 for the Earth and 10“ 6 for the 
•.Sun, so that solar system experiments are much more promising. 

As was pointed out originally by Nordtvedt 15 , a differential 

acceleration towards the Sun of the Earth and Moon would result if 

gravitational self-energy lacked passive mass. The resulting 

polarization of the lunar orbit about the Earth is detectable by 

lunar laser ranging using the retroref lectors left on the Moon 

during the Apollo program. Results to date indicate no anomaly, 

the accuracy in the measurement of the Eotvos ratio of gravitational 

16 

energy being about 2% 

If this effect existed, the Earth would exhibit an anomalous 

• 1 9 

acceleration towards the sun of order 10 g. While this could in 
principle be detected by a sensitive accelerometer in a laboratory, 
difficulties with first- and second-order solar gravity gradients, 
tidal effects, seismic noise, etc., probably preclude such an experi- 
ment. 
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It would be quite surprising if this type of Eotvos anomaly 
were discovered. One third of the anomalous precession of the 
perihelion of Mercury can be construed as arising from the active 
irtass of energy stored in the solar gravitational field^ 0 (there 
is more spherically symmetric mass inside the planet's position at 
aphelion than at perihelion) . Moreover , if the passive mass of a 
gravitating system changes when gravitational energy is converted 
to other forms (e.g., kinetic), it is clearly possible to conceive 
gedanken experiments which would exhibit small violations of energy 
conservation. 

A somewhat similar effect can arise in those theories which, 
while obeying the Equivalence Principle, allow spatial variations 
in the gravitational constant (e.g., the Brans-Dicke theory 3 ). 

Because the gravitational self-energy depends on G, spatial gradients 

© 

of G can produce small anomalous forces on a massive body-^. 

1.2.3 Spin-Orbit Interactions 

It is predicted by general relativity- 1 - 8 that a spinning body 

does not follow exactly a geodesic in the Riemannian space-time 

determined by neighboring bodies, when it is in acted on by no 

forces. In other words, a gyroscope in a gravitational field experiences 

• 19 

an anomalous acceleration, which has been calculated by Schiff 
for the case of a spherically symmetric, static field as 

a = 3GM f (r *H) (r x v) + (r*v) (r x H) ] [5] 

“ 3 2 5 

.. mc^r 3 
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where M is the mass of the source of the field, m is the mass 
of the gyro, r is the radius vector to the gyro, H is its angular 
momentum, and v is the orbital velocity. For a circular orbit, for 
which r • v = 0, with H lying in the plane of the orbit, an Eotvos 
experiment which compared the gyro with a non-spinning body would 
reveal an Eotvos ratio which varied sinusoidally at orbital period, 
of amplitude 
3Hft 


where ft is the orbital angular velocity. For example, 
a large gyro in low Earth orbit, with a wheel radius of one meter, 
spinning at 10,000 rpm (close to the value at which centrifugal 
stresses would tear it apart) , would exhibit an Eotvos ratio of order 

-17 

10 . Careful design and a sufficiently large apparatus might 

make this effect experimentally detectable, thereby providing 
another test of general relativity. Although such a device would, 
in principle, be capable of detecting the difference between an 
inertial acceleration and a gravitational field, this would not 
necessarily violate the Strong Principle of Equivalence, which is 
limited to infinitesimal regions. A spinning body cannot be in- 
finitesimal, even in principle, essentially because the periphery 
must move at less than the speed of light. Thus this experiment 
would constitute a test of the consequences, and not the foundations 
of general relativity. ; 
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1.2.4 The Electromagnetic Radiation Reaction 

At first sight, radiation from an accelerated charged particle 
seems to violate the Equivalence Principle. How does a charged 
body, at rest in a terrestrial laboratory, know that it is experiencing 
a gravitational field and not an acceleration, so that it does not 
radiate? The usual answer to this, of somewhat dubious validity, 
is that radiation phenomena must occur over distances of at least 
a wavelength, and hence cannot be considered as an experiment in an 
infinitesimal region. In any case, it is of interest to calculate 
the anomalous Eotvos ratio which would be expected if a charged 
body were compared with an uncharged one in Earth orbit, because of 
the radiation reaction. 

The power radiated by an accelerated charged particle is given 
by the Larmor formula^O 

P = - mxv 2 = £ rad *v [7] 

where 

2 ©i 

T = 2 me 2 ^8] 

where e is the charge. In a circular orbit, 

v 2 = (vxg) 2 = v 2 n 2 E9] 

so that [7 ] may be written 

!■— rad = ra ^ 2 v] ~ 0 [10] 

i ; —■ i ■ i —r 1 ' ? I < . - i 


I 


] 




or 

F = (R x fl) = - mg x xft 

-rad — ~ — 

where R is the radius of the orbit. The radiation 

reaction force is thus a drag, opposite to the orbital velocity, 

whose magnitude is less than the weight of the particle by the 

factor xfi. The characteristic time r has its maximum value, 

-24 

6.26 x 10 seconds, m the case of an electron, which, in low 
Earth orbit, would thus exhibit an Eotvos ratio of at most about 
10“ 26 . The effect is therefore negligible. 
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CHAPTER TWO 

AN EOTVOS EXPERIMENT IN ORBIT 


II. 1 Basic Considerations for a Space Experiment 

As discussed in the previous Chapter, there is now good 

| • . ' . 

evidence that all the forms of energy which make up matter obey 

the Equivalence Principle, with the possible exception of energy 

stored in the weak interaction. It appears that this last form 

of energy can be checked by carrying out an Eotvos experiment 

with a sensitivity improved by several orders of magnitude, to 

14 

about one part in 10 

The advent of space technology has made possible a substantial 
improvement in the experimental conditions for the Eotvos experiment. 
In low Earth orbit, there exists a balance between gravitational 
and centrifugal accelerations whose magnitude is close to one gee. 

.. .. c r 

This is to be compared with the forces available to Eotvos - 
(1.7 milligee) and Dicke et al 6,7 (0.6 milligee). Moreover, the 
free-fall environment in orbit allows great reduction in the problems 
associated with suspension of the apparatus, such as coupling to 
ambient noise. In principle, then, the Eotvos signal can be increased 
by three orders of magnitude, and the noise level substantially 
reduced. In view of the importance of the experiment to the 
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foundations of general relativity, it is an obvious candidate 
for performance in space. As discussed in Section III. 3, 
it may be an ideal experiment in terms of making use of the 
capabilities of the Space Transportation System (space shuttle) , 
in the development of the apparatus as well as in performance of 
the final experiment. 

The most obvious technique for carrying out this experiment 
in orbit is to use an adaptation of the rotational balances which 
have been successful in terrestrial experiments. However, gravity- 
gradient torques may swamp those due to Eotvos forces unless great 
care is taken to make the system inertially symmetrical. Without 
going into the details of the design of such a balance, some general 
conclusions may be drawn. 

For simplicity, consider a system with one of its ‘ principal 
axes along the orbit normal. The gravity-gradient torque is then 
along the orbit normal and of magnitude (calculated in Section III. 2. 4). 

T = - isi 2 Al sin 2 0 [12] 

g ^ 

where Al is the difference in the moments of inertia 
about the principal axes which lie in the orbit plane and 0 is 
the angle between one of these axes and the local vertical r 

Let us suppose the system is constructed of two different 
materials, A and B, for which the ratios of passive to inert 
mass are k (A) and k (B) , respectively. The Eotvos torque about 
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the center of mass is then 

= - £ x [k(A)/ A r dm + k(B),f B r dm] 

= - n£ x t/ B r dm] 

= - nm g x a [13] 




where £ is the local gravitational field, r is the 
radius vector from the CM of the system to an element of mass dm, 
and 


= r dm 

m B — 


[14] 


With appropriate symmetry, a will lie in the orbit plane, but 
generally not along one of the principal axes. is then along 

the orbit normal and of magnitude 


T e = - nmga sin 4> [15] 

where is the angle between a and the local vertical. 

In order to estimate the difficulty of inertially balancing the 
system, a reasonable condition to impose is that the frequency 
of gravity-gradient oscillations oe less than that due to Eotvos 
torques. If I is the moment of inertia about the axis along the 
orbit normal and b the radius of gyration about this axis, this 
condition may be written, using [12] and [15] , as 

M ^naR 

: I'.’ 3F 2 " [16] 
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where R is the radius 
geometry, the gravity-gradient 
size of the system than do the 
apparatus is required to allow 
a « b a 10 cm, [16] requires 

■— < 2 x lo" 7 


of the orbit. For a given system 
torques increase faster with the 
Eotvos torques , so relatively small 
adequate inertial balancing. For 


[17] 


-14 

in order to achieve an accuracy of ri * 10 . Inertial 

balancing to this accuracy is difficult but not impossible with 
current technology, especially as it may be possible to exploit 
the double-angle dependence of [12] and a design difference in the 
directions of u and a principal axis to effect final balancing 
on orbit. 


Even if this condition can be met, it must be recognized that 
the torques under consideration are extraordinarily weak. From 
[15], the angular frequency of oscillation due to the Eotvos 
torques is given by 




[18] 


which, under the conditions assumed in [17] , gives a 
period of some 2000 hours! It is clear that very great care 
indeed would be required to protect the balance from non-gravitational 
disturbing torques, such as those due to the Earth's magnetic 
field, radiation pressure or perhaps residual atmospheric drag. 
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Furthermore, it is highly desirable to design an experiment from 
which the data can be extracted more quickly. 


II. 2 Force-Balance Techniques 

In order to overcome the difficulty of inertially balancing an 
EGtvds torque-measurement apparatus for use in free fall, an alternative ““ 
technique is under study in the M.I.T. Measurement Systems Laboratory 21 ' 2 
(and also at the Hansen Laboratories of .Physics, Stanford University 23 ), 
in which the Eotvos forces are . to be measured directly, without con- 
verting them into torques. The design was motivated by the follow- 
ing considerations: 

i) Gravity-gradient forces may be minimized by placing the 
centers of mass of the bodies to be compared as nearly 
as possible at the same point in space. 

ii) The tensor properties of residual gravity-gradient forces 
may be used to distinguish them from the phenomenon 
under study. 

iii) The Eotvos force may be modulated by rotating the apparatus, 
allowing operation at higher frequencies and minimizing 
the duration of the experiment. 

iv) The design allows resonance to be used to enhance the 
sensitivity to Eotvos forces. 

v) Force-measuring devices (accelerometers) are simpler to 
instrument than torque-measuring devices, which must 
include gyros and/or star-trackers. 

vi) While a force-measurement experiment may seem, at first 

sight, more subject to disturbances due to external forces 
(magnetic, residual aerodynamic, etc.) than a torque 
balance might be, in fact the force balance technique 
allows discrimination against these effects in a way 
Which is difficult or impossible with a torque balance. 

In any case, if necessary these forces can be reduced 
or eliminated by well-established "pure gravity orbit" 
techniques. 24 ■/;/ 
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As sketched in Fig. II, the proposed apparatus is a satellite 
which consists basically of an aluminum wheel, spinning about its 
axis of symmetry, which is nominally aligned normal to the plane 
of the orbit. This orientation is, of course, stable under gravity- 
gradient torques. A sensitive electrostatic accelerometer is 
mounted radially in the plane of the wheel; it consists of a cylindrical 
electrode structure with an internal test mass constructed of 
material A and an external, annular test mass constructed of material 
B, as sketched in Fig. III. The accelerometer is so designed that, 
when the two test masses are at their null positions, their centers 
of mass are nominally coincident with each other and with the center 
of mass (CM) of the satellite. The object of the experiment is, 
of course, to compare the Eotvos ratios of materials A and B, by 
a differential measurement of the forces which they experience as 
the satellite moves along its orbital path, 

11.3 Analysis of the Experiment: Ideal Case 

Let m, M be the inert and passive masses of test mass A and 
m' , M' those of test mass B. The null positions are taken accurately 
coincident with the CM of the system and the displacements of the 
test masses from null are denoted by x, x' respectively. The mass 
of the wheel is assumed to be so large that motion of the test 
masses does not appreciably shift the position of the CM within 
the system. The satellite moves in a perfectly circular orbit 
about the Earth, which is taken to be spherically symmetric, and 




Figure III. Dual Accelerometer 
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no external forces other than gravitation act on the system. 

Ideally, the equation of motion of test mass A is then 

mx = ~mq+m w 2 x + [Mg_ -- mg_*] .i^ [19] 

where w is the inertial angular velocity of the wheel, 
q is a servo restraint force applied to the test mass, i_ is a unit 
vector along the (common) accelerometer axis, g is the gravitational 
field at the test mass and g* is the inertial acceleration of the 
CM of the system, due to terrestrial gravitation. In this equation, 
forces which are perpendicular to the sensitive axis, such as the 
Coriolis term, have been dropped: the cross-coupling between support 

and sensitive axes of the accelerometer is assumed zero. 

A similar equation to [19] , using appropriately primed quantities 
applies to test mass B. 

Without loss of generality, the Eotvos ratio of the satellite 
as a whole may be taken as unity. This is equivalent to defining 
the gravitational constant as that measured in a Cavendish experiment 
using masses of the same average composition as the system. Then 

g* = - a 2 R .■ [20] 

where is the orbital angular velocity of the satellite 
and R is the geocentric position vector of its CM. 

The gravitational field may be expanded in a Taylor series 
about the CM, most conveniently in tensor notation. For future 
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reference, we write the field for a general position r^ {i = 1,2,3) 
relative to the CM and compute the second-order terms, although 
only the first-order gravity gradients are required at present. 

For test mass A, 


Mg* - mg| = - g*] 


9 2 g, 


= m[Ck-l>gJ * + 2 T r j r k [505rl* + •••) 

3 J 3 k 

where k = M/m and * means the quantity is to be computed 
at the CM. Since k is very close to unity (as demonstrated by 
terrestrial Eotvos experiments) , it may be taken as such in [21] , 
except in the first term, so that 


[ 21 ] 


Hg i - mg i = Ck-l)g| * r..r. ♦ ♦ ] 


[ 22 ] 


r .here the first- and second-order gravity gradient tensors 
are given by the partial derivatives of the terrestrial gravitational 
field as 


r.. = 

9R T 




= - fi 2 [6. . - 3R - k.R.] 

1] A J 


and 


i. »*<i 


ijk 


2 L 3R. 3R, - J 
3 k 


* i# 2 t R i 6 jk + R j s ik + R k s ij - 


[23] 


[24] 
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where is the mass of the Earth and $ j_ j is the unit 

tensor. 

Substituting the first-order terms in [19] yields 

x + q - [to 2 - ft 2 + 3ft 2 R~ 2 (R.i x ) 2 ]x = (k-ljg*^ [25] 

or, singe i„ is rotating with respect to the local 
vertical (i.e., R) with an angular velocity (to-ft) , 

x + q - [to 2 + li2 2 + cos 2(to-fl)t]x = (k-l)g cos (to-fi) t [26] 

2 2 

If q is a simple spring restraint, 
q - Kx [27] 

this is recognized as a standard. Mathieu equation, driven 
by sinusoidal forcing function proportional to the difference 
between the EOtvos ratio for test mass A and that for the whole 
system. 

By writing out the equation similar to [26] for test mass B, 
using primed quantities, and subtracting from [26], we obtain 

y + Aq - [to 2 + i^ 2 + 1& 2 cos 2 (to-ft) t]y = ng cos (to-ft)t [28] 

2; 2 - 

where y = x-x 1 , Aq = -q', and n = k - k' is the EotVos 
number for materials A and B. A Mathieu equation for the differential 
motion of the two test masses results if q ! is also a spring restraint 
with the same spring constant per unit mass as q. * 
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The ideal performance equation 128] demonstrates the possibility 
of carrying out an Eotvos experiment with apparatus of this type. 

It would, of course, be possible to use an accelerometer with a 
single test mass, made of the material to be tested against the 
satellite as a whole, according to [26] . The advantages of the dual 
accelerometer are: 

i) Any pairs of materials may be compared, merely by changing 
proof masses. This flexibility may be used to enhance 
accuracy in the search for violations of equivalence by 
different forms of energy, as discussed in the previous 
chapter . 

ii) More importantly, as shown in the next Chapter, the con- 
figuration allows discrimination against many disturbing 
forces, which do not appear in the ideal equation [26]. 

It is to be noted that the amplitude of the Eotvos acceleration 
in [28] is Pg. in designing a force-balance Eotvos experiment with 
an accuracy of one part in 10^, a prerequisite is an accelerometer 
with a sensitivity (to sinusoidal acceleration) of 10”-^g. The 
problem of designing and, most importantly, testing such an accelero- 
meter is taken up in Chapter IV. 
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CHAPTER THREE 
SATELLITE SYSTEM DESIGN 

III.l Preliminary Stability Analysis 

In beginning a more detailed design of the force-balance 
satellite experiment discussed in the last Chapter, the first question 
to be taken up is that of determining the conditions under which 
Eq. [28] (or Eq.[26]) has stable solutions. To this end, we assume 
a servo restraint force between the test masses of the form 

Aq = Ky + 2 Ay [29] 

which represents a simple spring, with damping. The 
equation of differential motion [28] is then 

y + 2ly + [W 2 - cos 2-(m-fi)t]y = ng cos (w-ft)t [30] 

where 

W 2 = K - a ) 2 - isi z 
2 

For resonance, one might first choose X =0 and 
W 2 = (w-fl ) 2 


[ 31 ] 


30 , 


but unfortunately this ideal, infinite-Q case results in 
divergent solutions of the Mathieu equation 2 ^. The degree of damping 
required for stability may be calculated by the techniques described 
in Ref. 23. We consider the homogeneous version of [30] and, to 
put it in a standard form, first write 

t = (oj-fi)t [32] 

which yields (with ^ now denoted by y, etc.) 


y + y + — i [W 2 - cos 2x]y = 0 

(oi-n) (u-fi) 2 2 

The substitution 
y = e — Xt/(co — ^)z(t) 

reduces this equation to the standard form 
z + (a - 2b cos 2i)z = 0 
with 

a = (W 2 - A 2 )/ (w-0) 2 


[33] 


[34] 


[35] 


[36] 


and 

b = 2 [ 37 ] 

4 : . ■ . ;• . ■ ■ 

If [35] has stable solutions, then [33] certainly does also; 
but, because of the damping factor in [34], it is possible for [33] 
to be stable even though z(x) is divergent. We must therefore 
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investigate the properties of the unstable solutions of [35] , 
for different values of a & b. The stability diagram is Fig. 8 
of Ref. 25: the relevant region of the a-b plane, for operation 

: i- 

near resonance (a ^ 1) is shown in Fig. IV. The curves separating 
stable from unstable regions are given by 26 

a* = 1 + b ib 2 - -g-^b 3 - 1535 b 4 + ... 138a] 

and 

at = 1 - b - ^b 2 +^b 3 - “ ••• [38b] 

If the point (a,b) lies on the curve [38a] , the solution 
of [35] is periodic and even in t. If it lies on [38b], it is periodic 
and odd in x. Between these curves, the solution is divergent; 
it may be written in the form 2 ^ 

k=°° _ k=°° 

z(x) = Ae^ T l c k e 2 3^ T + Be“^ X I c^e -2 ^ 1 [39] 

k=-oo k=-°° 

where the summations represent periodic functions , A & B 
are constants determined by the initial conditions, and y is 
a function of a & b. It is clear from [34] that the solutions to 
the original equation [30] will be stable if 

X > (a j-fi) 1 ia 1 = X m±n [40] 

Between the curves [38a] and [38b] , we may choose a parameter 
tf such that 2 ' /'V‘ S ' : ■ 





[411 


a = 1 - b cos 2a + ib 2 (icos 4a - 1) + ^b 3 cos 2a 

4 2 64 

+ t ^ 4 <1 - i 002 4a) - • • • 

which reduces to [38b] for a = 0 ana to [38a] for a = tt/ 2. 
Given values of a & b, [41] may be solved numerically for a and 
then 2 ^ 

U = - ib sin 2a + -T-I^b 3 sin 2a - ^r- .. 3 .b 4 sin 4a - ... [42] 

Z ±Zo J.U24 

We consider two modes of operation of the system which might 
be used in practice: 

i) w>>n 

If the wheel is spinning at an angular velocity much larger 
than orbital angular velocity, [37] shows b<<l. It is then sufficient 
to retain only up to the quadratic terms in b in [41] . To estimate 
the required damping, we assume operation at resonance (a ” 1) , 
and then [41] becomes 

ib 2 (i- cos 4a - 1) - b cos 2a = 0 

or 

cos 2 2a - ^cos 2a - * 0 [43] 

to second order in b. 


which may be solved to yield. 
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To this accuracy, J42J then gives 


U » - 2b 


[45] 


.and,, from [40], the minimum damping for stability is 


ft' 


Wn “ 2 (W “ fi)b 8 <ui~ft) 


[46] 


Under these conditions, the time-varying coefficient in [30] 
produces a small perturbation to a simple second-order harmonic 
system, which, according to [46], has 


Q, 


(co-ft) _ 1 


max 


2A 


mm 


2 U 


_ 1. — j4 /to t \ 2 

b 3 v ft ' 


[47] 


It is thv-g clear that very light damping may be used if oo is 
high enough. For example, in low orbit, with co = 10 rpm, a value 


Q 


max 


10^ is obtained. 


If the solution [39] is substituted in [34] and t replaced by 
t, by [32] , it is easily seen that the effective damping factor is 


^eff ~ ^ ~ (w-ft) | ]i 


[48] 


so that 


^ _ (co-ft) 

Qeff “ 2A 

^ A eff 

can be much higher than C> „ if A is close to A • . 

lLlaX mill 

Comparison with an accurate numerical calculation shows that 


[49] 


[46] is an adequate expression for A 


mm 


if co is greater than about 


3ft 
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ii) a) = ,0 CInertially stationary wheel) . 

In this case, from [37] , b = 3.. If we put A = A m ^n and use 
the resonant condition [31] in [36], we find, using [40], 

a = 1 - u 2 150] 

Since we cannot be assured a priori that y is small, this 
expression must be used in [41] and the equation solved simultaneously 
with [42] . An iterative numerical procedure (reguli falsi) gives 

A = .365 [51] 

Note that the expression [46] gives a result which is accurate 
within about 3% even in this extreme case. While high effective 
Q can still be achieved in principle, in the presence of the 
heavy damping implied by [51], in such a case small changes in the 
damping could lead to instabilities. This is one of the reasons 
the spinning system is preferred for this experiment. 

In practice, of course, a simple spring with damping is a 
very elementary servo. A more realistic design requires study of 
disturbing forces acting on the system, and the question of system 
stability must be examined again in the light of this analysis. 

For example, it would be possible to compute and apply a correction 
for the time-varying coefficient in [30] , if this should prove 
desirable. However, [46] provides a good rule of thumb for estimat- 
ing the minimum allowable damping in the absence of such compensa- 
tion. 


T 
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III. 2 Disturbing Forces 

In the system under consideration the Eotvos acceleration, * 

if it exists, always lies along the local vertical, and is modulated 
by rotation of the wheel. Any forces applied to the wheel or to 
the tost masses which appear, in a reference frame attached to the 
wheel, at frequencies close to (w-fi) , the angular velocity of the • 

wheel with respect to a local vertical reference frame, may 
masquerade as Eotvos forces, although considerable discrimination 
is possible when the differential motion of the two test masses 
is regarded as the output of the system. Limits to this filtering 
action of the dual accelerometer will be considered below, but, * 

to maximize sensitivity, it is clearly desirable to minimize extraneous 

ar 

forces, in particular those, such as residual atmospheric drag, 
which are approximately constant in the orbital frame. Those 
forces which have constant components along the local vertical may 
prove especially troublesome. » 

Techniques exist^ 8 for isolating the system from most external ; 

forces, by making it "drag-free". Indeed, it was the availability 
of such techniques which provided one motivation for designing 
the system to measure Eotvos forces directly. However, only those 
forces which, directly or indirectly, appear. along the sensitive 
axis of. the accelerometer are of interest, and this may allow a 
simpler mechanization than that involved in forcing the entire * 

satellite to follow a pure gravity orbit. This possibility of 

* 

course depends on the feasibility of designing an accelerometer 
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with sufficiently low cross-coupling between the support and 
sensitive axes. 


The list of possible disturbances which must be taken into 
account includes at least the following effects: 

1. The orbit of the satellite is slightly elliptical. 

2. The center of mass of the system does not coincide with 
the center of gravity. 

3. The system moves in the vicinity of a massive object, 
which exerts gravity-gradient forces on it (for example, 
the experiment is free-flying in the shuttle payload bay) . 

4. The wheel angular momentum is not along a principal 
axis, nor along the normal to the orbit plane. 

5. The null positions of the test masses do not coincide 
with each other, nor with the CM of the system. 

6. The sensitive axes of the dual accelerometer are not 
coincident, nor are they perpendicular to the spin axis 
of the wheel . 

7. Support forces for the test masses couple into the sensitive 
axes, due to Coriolis effects as well as to mechanical and 
electromagnetic imperfections. 

8. The CM of the system shifts when the test masses move. 

9. The test masses are affected by stray electromagnetic 
forces , generated by the satellite or arising from natural 
ambient fields. 

10. Various mechanical forces (noise) are applied to the 
wheel. 

11. The displacement detector is noisy, exerts forces on the 
test masses, and possibly exhibits non-linearities 
(e.g., threshold). 

System design for this experiment is not yet complete, and 
not all of these effects have as yet been taken into account. 


111.2*3. Orbit Eccentricity 


If the orbit is not precisely circular, the first-order gravity 
gradient tensor [23] must be written 

r ij = - GM e R~ 3 [S i;j - [52] 

where jj_ is a unit vector along the local vertical, 
to allow for variations in the orbit radius R. The equation of 
motion [25] of the first test mass is then 

x + q - [w 2 - G%R“ 3 (1 - 3(i.ix) 2 ]x = (k-l)GM R" 2 j.i x [53] 

If we assume for simplicity that the accelerometer sensitive 
axis is along the local vertical at perigee, then 

i.i^ = cos (tot - (j)) [54] 

where cf) is the true anomaly, the instantaneous geocentric 
angle away from perigee in the orbit. The radius of the orbit is 
given by 2 9 

R - E r [55] 

1 + e cos <j) 

where e is the eccentricity and p is the semi-latus 
rectum, which is the radius of the orbit when the true anomaly 
is 90°. 

In order to find the variation of (j) with time , we first 
introduce a new variable E, called the eccentric anomaly, by 
the relation 


The eccentric anomaly is the solution of Kepler's equation-^ 

E - e sin E = 0 , o t [57] 

where is the mean orbital angular velocity. This 
equation may be solved in terms of a Fourier-Bessel series 
expansion for E, but, to first order in e, the solution is obviously 

E - ft Q t + e sin fi Q t [58] 

If we write 

(ft = a o t +5 [59] 

and insert this and [58] in [56] , we find, after 
some reduction, to first order in Sand e, 

5 = 2e sin ft Q t [60] 

and then [54] gives 

= cos(w-Q 0 )t + 2e sin ft Q t sin(w-fi 0 )t [61] 

Now, to first order, 

GM 0 R" 3 (1 - Sti-i^) 2 ) = fl2(i + 3e cos ft Q t) [1 - 3 cos 2 ( u -^)t 

-* 12e cOs(w“^o) t sin ft Q t sin .Cu>“fl 0 ) t 
= Q2[- i - | cos 2 ( 0 J-S 2 o )t 

+ 3e[- — cos ^ 0 t - 2- cos (2 co - 3^,)t + 2. cos (2w-fi 0 ) t] ] [62] 
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and, on the right side of [ 5 3 j , 

= g Q (l + 2e cos ^ 0 t) (cos ( w ~' fi o) t + 2e s ^- n fi ot sin(03-fi o )t 
= g Q [cos (to-fi 0 ) t + 2e cos (u)-2ft Q ) t] [63] 

Suppose now that there is a steady state offset Ay Q between 
the null positions of the two test masses. When the equation of 
differential motion is formed, analogously to [28] , by subtracting 
from [53] a similar equation for the other test mass, the first 
order effects of the orbit eccentricity will be, firstly, to 
change slightly the time-varying coefficient, according to [62] ; 
and, secondly and more importantly, to introduce a driving force 
which, apart from constant terms, is given by 

f = fi 2 Ay 0 [A cos 2(o)-^ 0 )t + 3e (i cos fi Q t + Z. cos (2oo-3fi 0 ) t 

- cos (2aj-fi 0 ) t) ] [64] 

4 

Note that there is no component of this force with the principal 
frequency (w-fi) at which the Eotvos acceleration appears. In 
fact, the double harmonic term (independent of the eccentricity) 
can be used to drive the null positions of the two test masses 
into coincidence. 

If the wheel is inertially stationary, however, a term 

f l = Ay 0 C0S : ft 0 t ; : A; An [65] ■■■: 
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appears in [64] , at the same frequency S2 Q exhibited ' 
by the Sotvcs acceleration in this case. Since the amplitude of 
the Eotvds signal is ng = nft^R 0 , we must have 

Ay 0 << 4 . n . k o [66] 

3 e 

if f-j_ is not to be falsely interpreted. For e - .002, 
a reasonable value for low Earth orbit, and if the design sensitivity 
is ii = 10“ the null offset must be small compared to 40 microns. 
This is not a difficult requirement, but the effect could become 
troublesome if it were desired to attempt much greater sensitivity. 

In any case, the problem can be avoided by operation in the spinning 
mode. 


III. 2. 2 CM-CG Deviations 


Consider an arbitrary rigid body in circular orbit. Set up 
a reference frame with origin at the center of the Earth, rotating 

' . ; j f - 

at orbital angular velocity so as to keep the 1-axis along the 
orbit normal and the 3-axis along the local vertical. In this 
frame, the total (gravitational plus centrifugal) force acting on 
an element of mass dm in the body, located at a position vector 
with respect to the CM, is (cf. Eqs. [23] and [24]). 


dFj_ = - fi 2 (% + rjjdm + g*dm + T^jrjdm + T^j^rjr^dm + . . . [67] 


where is the geocentric position vector of the CM 
and g* is now the gravitational field at the CM. By definition 
of the CM, 
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/r^dm =0 [68] 

so that integration of [67] over the body yields 
F i = m(g£ - + T i j k fr^rydm + ... [69] 

The definition of the inertia tensor of the body is 2 ® 

!jk = ^ (r n r n«jk ' r j r k )dm [701 

so that the last term in [69] may be written 

AFj_ - ^ r n r n^ m ~ ^ijk^-jk 

=: ■” ^ijk^jk 

= “ + 211.1^ - 5R" 2 R i R j R k I jk ] [71] 

since, from [24] 

T ijk S 'jk = T ijj = §ii 2 R _2 [3Ri + Ri + Ri - 5%] = 0 [72] 

Now let us suppose that the principal 1-axis of the body lies along 
the orbit normal, and that the principal 3-axis makes an angle 0 with 
the local vertical. In the principal-axis frame, Ij^ is of course 
diagonal, and 

R^ = R[0, sin 0, cos 0] [73] 

Rjl i> = R [0 , I 2 sin 0, I 3 cos 0] VV; -h [74] 
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and 


RjR^I j k = (I 2 sin 2 e + I 3 cos 2 0)R 2 


[75] 


so that [71] may be written, after seme trigonometric 
reduction , as 


AF - 


■•?mgR“ 2 

2 


<*! 


£ £ 2 

4 ^2 

I 

4 ^2 


i £ 2 ) sin 0 (£ 2 -£ 2 )sin 3 e 

^ 2 ) cos 0 +|. (? 2 -? 2 )cos 30 


[76] 


where , £ 2 ' £3 are the radii of gyration. Compared 'to 
an infinitesimal body, a finite-sized body in orbit thus, in 
general, experiences a small additional gravitational acceleration. 
We note first that, if 0 = 0, the incremental acceleration is 
along the local vertical, and of magnitude 


gR" 2 df + - 2£j) 


[77] 


If 


£ 2 = £ 2 as in the case of a wheel spinning about the normal 


to the orbit, then, as is to be expected, [76] shows that the accelera- 
tion is along the local vertical, regardless of 0 4 If the wheel 
is taken as a uniform disc of radius a, then [ 77 ] gives 




For low Earth orbit (R - 6550 km) and a si 50 cm. 


[78] 


e - 2 x 10' 


•15, 


[79] 


The test masses are physically much smaller, so the acceleration 
applied to them is expected to be at least two orders of magnitude 
smaller than [79] , depending on their design. Thus [78] appears 
in the system exactly like an Eotvos acceleration, as if the ratio 
of passive to inertial mass of the wheel were greater than unity, 
except that it is sensed by both test masses. The filtering 
capability thus provided should permit, as far as this effect is 
concerned, operation down to about q = 10 - -'-’ 7 . At this level, it 
becomes necessary to take into account the differences in the 
detailed inertial properties of the two test masses -- specifically 
the difference obtained when [76] is applied to each of them. 

The effect [79] could perhaps be reduced by an order of 
magnitude by making the entire apparatus as small as possible. 

A better technique would be to make the "wheel" isoinertial, as 
[76] and [77] vanish if the radii of gyration are all equal. In 
this case, however, the system would not be positively stable in 
its nominal attitude under gravity-gradient torques, and a separate 
attitude control system would be required. 

Of course, it is not necessary that the system be isoinertial 
in order for [77] to vanish. If we choose 

Cl = 2^3 - Cl [80] 

then it is easy to see that the principal 1 -axis will have 
maximum moment of inertia, as desired for attitude stability, as 
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long as These conditions are met, for example, by a flat 

rectangular plate, lying in the plane of the orbit with its long 
dimension horizontal, if the length-to-width ratio is y2:l. However, 
the fact that the acceleration vanishes for particular orientations 
of particular asymmetric bodies is academic, since it is necessary 
to rotate the apparatus* to modulate the EQtvds forces. What is 
significant is the component of [76] along the sensitive axis of 
the accelerometer, with 0 = (ou-fi)t. If, for example, the accelerometer 
is along the principal 2-axis, it would be desirable to have 



3 

4 



[81] 


since the third-harmonic term may be filtered out. From 
this, however, the condition implies > £ 2 *, whereas 

> K 2 implies Thus [81] is incompatible with being 
the maximum radius of gyration; the best that can be achieved is 
inertial symmetry. 


III. 2.3 Nearby Massive Bodies 

For a preliminary evaluation of the feasibility of operating 
the experiment in the vicinity of other equipment {e .g. , free- 
flying in the payload bay of the space shuttle), we consider here' 
the gravitational effects on the system of a body of mass M° , located 
at a vector position R° with respect to the CM of the system. For 
numerical estimates, we take M° — 1 ton, R° = 2 meters. 

♦The apparatus is, of course, rotating relative to the orbital frame 
even in the inertially stationary case. 


At the CM of the apparatus , the gravitational field of this 
object is of order 10”^g. The presence of the body thus makes a 
negligible change in the Eotvos signal (the right-hand side of 
[28]). The direct effect of this field, the acceleration of the 
apparatus towards the body, is also negligible, being two orders 
of magnitude smaller than the apparent relative accelerations due 
to the terrestrial gravity gradient. 

The first-order gravity gradients due to the disturbing body 
can be more significant, changing the time-varying coefficient in 
[28]. Close to the surface of a spherical body of density p, the 

># 

gravity gradient tensor will have components of order 

GMR" 3 k i ttGp [82] 

3 

and thus can be comparable with the terrestrial gravity 
gradients if the density is comparable to that of the Earth. We 
can in fact write [28] as 

y 4 - Aq — [to 2 + (1 + 3 cos 2(to-fi)t) 

2 ® 

+ i.GM°R°‘~ 3 (l + 3 cos (tot -0))]y = ng cos (w-f2) [83] 

2 

where 0 is the angle between R° and an inertial reference direc- 
tion and it has been assumed for simplicity that R° lies in the 
plane of the orbit. If R°.R is constant (e.g. , if the shuttle is 
stabilized with respect to the local vertical) , the time dependence 
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of the disturbing gravity gradients will be the same (apart from 
a phase angle) as that of the terrestrial term. In this case, the 
principal effect is an increase in the coefficient b in [35] ; by 
[47], this may increase somewhat the damping required for stability. 

With the numerical values for M° and R° given above, the disturbing 
gravity gradients are two orders of magnitude smaller than the 
terrestrial term. The conclusion is thus that, if it is desired 
to operate the experiment in the shuttle, a more careful study of 
the effects of first-order vehicle gravity gradients is warranted, 
but it appears probable that the sensitivity of the experiment will 
. not be compromised by this effect. 

Further difficulties arise when the second-order gradients 
are considered. If R° is in the plane of the wheel, whose radius 
is a, then the acceleration of the wheel towards M° , due to the 
difference in position of the CG and the CM, is given by [78] as 

£ -- ~ 4 x io”-*-lg [84] 

8 

with the numerical values and a = 50 cm. Comparison with 
[79] shows the effect is three orders of magnitude larger than 
that due to the second-order terrestrial gradients. 

The acceleration [84] is, of course, extraordinarily small 
by ordinary standards. If the experimental apparatus has a mass 
of 100 kg, for example, the force exerted on it due to this effect 


is only 4 millidynes; this is equal to the radiation pressure 
exerted by reflection of a light beam whose power is 6 watts! 

It is not possible to avoid the effect [84] by attaching the 
experimental apparatus in a suitable bearing to the structure of 
the vehicle (even if the mechanical noise problems thereby introduced 
could be overcome) , because the test masses would then sense the 
direct gravitational field of the disturbing mass (M.0“9g) as well 
as terrestrial gravity gradients ('vlO”’ 7 g if the attachment point 
is 2 meters from the CM of the shuttle). Free flight is essential. 

The acceleration e is similar to an Eotvos acceleration but of 
course does not represent a limit to the accuracy of the experiment 
if carried out in a large vehicle in orbit. In the first place, 
the effect is sensed equally by both test masses; the discrimination 
afforded by the differential measurement depends on the design of 
the dual accelerometer, but is expected to be several orders of 
magnitude. Further discrimination is possible if the vector R° 
varies in direction or magnitude, with frequency components different 
to those involved in the Eotvos measurement. For example, a slow 
rotation of the shuttle would remove the acceleration e from the 
passband of the accelerometer. In order to avoid the essentially 
continuous thrusting which this might involve if the experiment were 
not quite close to the CM of the vehicle, a possible technique would 
be to align the roll axis of the shuttle approximately along the 
local vertical and allow it to execute dumbbell gravity gradient 


oscillations, of an amplitude sufficiently small so that the apparatus 
does not strike the walls of the payload bay. In this case, there 
would be a constant component of e, but measurement of the component 
at the dumbbell frequency (in the sum motion of the test masses ) 
would allow computation and subtraction of it in data processing. 

More Herculean measures are possible to counteract this effect. 
For example, measurement of the acceleration, as discussed above, 
would provide the data needed to exert a compensating force on the 
apparatus (perhaps by radiation pressure!). Alternatively, it is 
in principle possible to add compensating masses, attached to the 
shuttle, around the apparatus, in such a way as to reduce the 
second-order gravity gradients . 

Without going to such lengths, it is probably possible to carry 
out the experiment in the shuttle payload bay, with an accuracy of 
perhaps one part in lO 1 ^*. If higher accuracy is required, launching 
the apparatus as an independent satellite is probably the preferred 
technique, paying the price for the increased autonomy which would 
then be needed . 


*At this level, It may be necessary to take into account second- 
order gravity-gradient forces on the test masses themselves. 
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III. 2. 4 Attitude Motions of the Apparatus 
a) Torque-free Nutation 

To start with the simplest possible case, we assume that the 
system is inertially symmetric about the spin axis, having a moment 
of inertia I about that axis and 1^. about transverse axes, and 
that the test masses are so small that their motion does not affect 
the motion of the system as a whol£i. The motion is then that of 
a classical rigid body. For future reference, we start by sketching 
the well-known derivation of the nutation frequency, in the absence 
of gravity-gradient or other torques^. 

The equation of motion is determined by conservation of 
angular raomen turn : 

[H] = [H] + oj x H = 0 [85] 

- I - B ~ “ 

where [H] is the rate of change of angular momentum relative to 
1 inertial space; 

[H]g is the rate of change of angular momentum relative to 
body axes; and 

ox is the inertial angular velocity. 

In body axes, 

w = [ax., ax , ax 1 , . [861 

— s 2 3 


= [IgU s , I t UY 2 , I t ax 3 ] 


and [85] reads, in components. 


w s = 0 

+ u n w 3 = ® 
u) 3 “ 4J n w 2 = 0 
where 

W n — [is ■” 1] Ug 
X t 

Differentiating [ 88 c] and substituting from [ 88 b] yields 
. 003 + = 0 

with the solution 
co 3 = A cos to n t 

which, when used in [ 88 b] , yields 
co 2 = ” A sin co n t 


[ 88 a] 

[ 88 b) 

[ 88 c] 

[89] 

[90] 

[91] 

[92] 


Let us now suppose that the 2-axis is along the accelerometer 
axis. Due to the nutation, the test mass experiences a centrifugal 
acceleration 

a c = - W 3 X = - A 2 x cos^ u) n t = - i. A 2 (1 + cos 2co n t)x [93 ] 

If the accelerometer is a simple spring with damping, its equation 
of motion (neglecting gravity gradients) is then 


1 


L I 
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x + 2 Ax + [ (K - o)| - ~ A 2 ) - i A 2 cos 2w n t]x = ng cos(w-fi)t [94] 

so that, in many respects, the effects of nutation on 
the system are similar to those of gravity gradients. If the 
amplitude of the nutation is variable, however, it shifts the 
resonant frequency of the system, which limits the allowable Q; 
if the nutation is uncontrolled, we must have, in fact, 

A < (w-Sl) (2/Q eff ) 1/2 [95] 

in order to stay within the 3-db bandwidth of the accelerom- 
eter, assumed resonant at (w-fi) . For reasonable values of Q e ff , 
this condition can be achieved during start-up of the experiment 
by the use of nutation damping devices. 

The stability of the test masses under nutation may be in- 
vestigated by the methods of Sec. Ill . 1 , using the homogeneous 
form of [94] . The equation can be reduced to the standard form 
[35] , where now 

a = (W 2 - A 2 )/o)2 [96] 

with 

W 2 = K - a) 2 - J. A 2 [97] : 

s 2 

and 

b '[-£-) 2 [98] • 

2o n 


Assuming the damping is small, the system is tuned to the 


Edtvos signal. 


a ~ (w-ft) 2 /w 2 


n 


[99] 


and, from [95], 


b < 2 /q 

2 w /w eff 
n 


2Q 


[ 100 ] 


eff 


According to [89] , choice of the ratio of moments of inertia 
allows control of the nutation frequency a) n which appears in the 
time-varying coefficient in [94], As w n is varied, according to 
[100] the operating point in the a-b plane remains closer to the 
a-axis than a line through the origin, of slope For any 

reasonable Q e ff/ examination of the stability diagram from the 
Mathieu equation 2 -* shows that the solutions of [94] will be stable 
except for quite restricted regions around 


a = n 2 


[ 101 ] 


where n is any integer. In other words, the condition 
[95] is sufficient to ensure stability as long as 

<o n f i(u)-n) [102] 

n 

or, from [89] (since io s » w) 

I /I. / 1 + ~(1 - -) [103] 

s' t . n ur 

where n is close to integral. 


The result [102] merely states that the nutation frequency 
should not be close to a submultiple of the resonant frequency of 
the accelerometer. This imposes some constraints on the design 
of the apparatus: for example, the wheel should not be closely 

similar to a thin disk in its inertial properties if w >> ft, 
for then I /I. = 2 and approximate equality will hold in [103] for 

o U 

n = 1 . 

b) Gravity-Gradient. Torques 

The first-order terrestrial gravity gradients exert a torque 
on the satellite which may be calculated as follows: In an orbital 

inertial frame (oi-frame) , which is an inertially non-rotating 
frame whose origin is comoving with the CM of the satellite, the 
net force on an element of mass dm of the satellite, located at 
vector position r^ with respect to the CM, is (cf.[23]) 

dFi = T^jrjdm = - ft 2 [r^ - 3R~^Rj_r.j_Rj ] dm [104] 

where it is assumed the satellite is in a circular orbit. 
In vector notation, 

dF = - ft 2 [r - 3R~ 2 (R.r)R]dm [104a] 

The total torque exerted on the body by the gravity gradient 
is the integral over the body of the moment of the force [104] 
(unless otherwise noted, all integrals below are over the body): 
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M = /r x dF 

= 3ft 2 R -2 / (R.r) (r x R)dm 

= - 3fl 2 R" 2 R x J [105] 

where we have used the fact that r x r = 0 and where 

J = / (R.r) r dm [106] 

Reverting now to tensor notation, we use the definition 
[70] of the inertia tensor to write 



/r^r^dm = Rj [/r 2 Sjj_dm 




[107] 


The first term on the right, in vector notation, is R/r^dm, 
which will vanish when the cross product with R is taken in [105]. 
It is therefore sufficient to consider the vector 


J i - " Vii 


[108] 


In principal axes i_, 2 , k in the body, this vector is 


1-jR.i, i - I 2 R ? i 2 " 1 IS 


[109] 


so that [105] gives 


M = 3n 2 R^ 2 [I ;L (R.i) (R x i ) + I^R^ (R x 2 ) + ^(R.k) (R x k) ] [110] 
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Assume now that the body has an axis of symmetry, so that 



and then 


* 

[ 111 ] 


M = 3fi 2 R“ 2 [I s (R.i) (R x i) + I t R x ((R.j)j_ + (R.k)k)] 

= 3ft 2 R“ 2 [x s (R.i) (R x i) + I t R x (R - (R.i)i)] 

= 3fi 2 (I g - I t )R" 2 (R.i) (R x i) [112] 

The gravity-gradient torque is thus always perpendicular 
to the plane containing the symmetry axis (i_-axis) and the local 
vertical, and its magnitude is 

M = |fi 2 (I s -I t )sin 20 [113] 

where 0 is the angle between the symmetry axis and the 
local vertical. 


Average torque ; 

The spin angular momentum is I 0 w_; comparison with [113] shows 

S b 

that the precessional angular velocity w will be smaller than 
the orbital angular velocity, approximately in the ratio of to 
co s . In this case, components of the torque at orbital and twice- 
orbital frequency are largely filtered out by the gyroscopic 
dynamics. While the effects of these higher-frequency components 





should be examined in a more complete analysis, because of the 
possibility that they might excite resonances which could be 
mistaken for Eotvos forces, attention is restricted here to the 
effects of the torque averaged over an orbit. 

To compute the average torque, set up a triad of inertially 
fixed unit vectors l_,m,n, with n along the orbit normal. Let the 

I 

components of the symmetry-axis unit vector JL along these directions 
be i^,i 2 /ig: these will be assumed not to change significantly during 

a single orbit. Measuring the geocentric angle <j> from the instant 
when 1^ is along the local vertical, we calculate the integral 

^ f Q (Xr.x) ( ia x 

where 


ip = 1 cos tp + m sin <f> 


[114] 


is a unit vector along the local vertical. Now 


i^.i^ = i^cos <J> + ^sin <p 


[115] 


iR x i = 


m n 


cos <j> sin <jj 0 


x 2 x 3 


i^sih 4 1_ - i^cbs <j> m + (^cos <j) - i^sin <j>)n [116] 


i 
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The integral is then 

1 2TT 2 

f o {(i^i 3 sm <j> cos <f> + i 2 i 3 sxn <f>)l 

2tt 

- (l^i^cos^cj) + (J) cos <j))ra 

2 2 2 2 

+ [.!•. i 9 (cos^tj) - sin^) + (i - i ) sin <j> cost|)]n} d<f> 
which reduces readily to 

i.[i 2 i 3 l - i'ii 3 £^ = ^i 3 (i x n) = |.(i.n) (i x n) [118] 

all other terms integrating to zero. 

This integral can now be used to find the average value of 
the torque [112] over an orbit: 

M = (I - I.) (i.n) (i x n) [1191 

— 2 S t — — ; ; 

The average torque is thus perpendicular to n (i.e., it lies 
in the orbital plane) and is of magnitude 

M = hi 2 (1-1. ) sin 2y [120] 

4 st 

where y is the angle between the symmetry axis and the 
orbit normal. 

c) Gravity-Gradient Induced Precession and Nutation 

The relationship between the set of unit vectors i_, along 
the principal axes of the satellite and the set l,m,n which are 
fixed in the oi-frame can be described by three Euler angles, 
as shown in Fig. V: 


[117] 







JL 
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Y , the angle between the symmetry axis i_ and the orbit 
normal n; 

a, the azimuth angle of the symmetry axis about n; and 

3 , the rotation angle about :L . 

The order of rotations from the oi-frame to the body frame 
(b-frame) is first, a rotation a about n, a rotation Y about the 
displaced 2 " ax i- s ' and finally a rotation 3 about the displaced 
i-axis . 

In body axes, the angular velocity of the body is given in 
terms of the rates of change of the Euler angles by 


w 


b _ 


• 3 

3 + a cos Y | 

• • . I 

Y sin 3 - a sin Y cos 3. 

! 

* * 

Y cos 3 + a sin Y sin 3, 


0 ). 


iii. 


!u. 


[ 121 ] 


and the angular momentum is given by 


Hb 


*t w 2 

^“3 


[ 122 ] 


There is no component on the average torque [119] along i 
or n. The components of angular momentum in these directions 
are therefore conserved. In particular, from [121]* 


*a and b are not to be confused with the coefficients in the standard 
form of the Mathieu equations . 
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/ 


I U), 
si 


• ® 


I g (3 + a cos y) 


I b = constant 
s 


[123] 


In body coordinates. 


b 


n 


cos y 

« cos 3 sin y 
sin 3 sin y 


[124] 


so that the component of angular momentum along n may be 


written 


w 2 

H.n = I g b cos y + a I sin y = 


I s a = constant 


[125] 


from [121]. Thus 


a - [a - b cos y] 

sin y 


[126] 


where 


0 * V ! t 


[127] 


and 


3 = b a cos Y = b - | ^ [a - b cos Y] 

sin y 


[ 128 ] 


I i I 1 1 £ l 
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We have thus found that a and 3 depend only on the value of 
y. To determine the values of the constants a and b, the variation 
of y with time , we need another equation, which can be obtained 
from the conservation of energy. The rotational kinetic energy 
is 


K = |[I S <4 * I t (^ + U 2)] 

= j I s b I t (Y 2 + i 2 sin 2 y) 


[129] 


from [121] . The potential energy in the gravity gradient 
force field may be obtained by integration of the magnitude [120] 
of the average torque : 


V = /Im dy 

O 

= - cos 2y 


[130] 


The total energy is then 
E = K + V = constant 

An equivalent constant of the motion is 

d = ~ [E - ~I b 2 ] = y 2 + a 2 sin 2 y - \ (a-l)cos 2y 

1 . / s 4 


[131] 


[132] 


since b is constant. Substituting for a from [126] yields 


. 2 
sin y 


(a - b cos y) + x f2 2 (d-l) cos 2y + d 
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or 


2*2 2 
sin y Y = - o (a 

If we put 
u = cos y 


- b cos y) + 4 -^ 2 ( 0 -l) cos 2y sin y + d sin Y [^33] 


[134] 


the equation becomes 


«s 2 = f(u) 

- c(l - u 2 )(2u 2 - 1) ♦ d(l - u 2 ) - a 2 (a - bu) 2 

where 

c = | n 2 (a-l) 


[135] 


[136] 


The problem has thus been reduced to quadratures , in terms of 
elliptic integrals. Some useful information may however be derived 
without carrying through the general integration. For example, 
the amplitude of the motion in y ma Y t> e found by noting that, at 
the limits of such motion, u is zero, and therefore the limits are 
given by the roots of f(u). 

To simplify the discussion, we consider the motion which 
develops if the spin axis is initially inertially stationary? 
i.e., we choose the initial conditions a == 0, y - 0 at t = 0 . 

Let the initial value of u be u Q . Then, from [126] at t = 0, 

a = bu Q [137] 





and then, since u is zero at t = 0, [135] gives 
d = - c(2u 2 - 1) 

so that, at general times,[135] may be written 

u 2 = f(u> = 2 c(1-u 2 )(u 2 -u 2 ) - 0 2 b 2 (u-u ) 2 

o o 

= Cu-u o )[2c(l-u 2 ) (u+u o ) - or 2 b 2 Cu-u o )] 

Apart from u Q , the roots of f (u) are those of 
2c(l-u 2 ) (u+u Q ) - 0 2 b 2 (u-u o ) = 0 

With the above initial conditions 
b = w 

s 

the spin angular velocity, so that 

2c/o 2 b‘ i = j— 2 . (0-1) /a 2 « 1 

s 


[138] 


[139] 


[140] 


[141] 



if the spin angular velocity is much greater than the 
orbital angular velocity. An approximate solution to the cubic 
[140] may be found by writing it in the form 

. u " u o = ^72 (1 " u2 Ku + u ) [142] 

0 b 

According to [141] , the quantity on the right is small, so 
that the root is close to u^, which means in turn that u may be 
replaced by u Q on the right, to give the approximate root 


U 1 * u ofl * 


[143] 
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The other two roots of [140] are imaginary, 


If now we write 


u = cos Y 
o 'o 


U 1 " cos (Y 0 + «j) a cos y o - 6 lS in y c 


[144] 


and insert these into [143] , we obtain 


f 2y o 


[145] 


for the amplitude of the gravity-gradient induced nutation 

(i.e„, fluctuation in y) • The effect is very small indeed, less 

( 

than 0.1 seconds of arc if u>_ = 10 rpm, in low orbit. 


The nutation frequency may be obtained by writing 


u = cos (y +6) = cos y - 6 sin y 


[146] 


When this is inserted in [139] , we obtain, to second order 


in 5 , 


( 6 2 = - 6[A 1 +■ A 2 S] 


[147] 


• ?■? 

i f-S 


where 


A. '= 2 c sin 2 y 


[148] 
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A 2 = ° 2b2 + 2c C 5u o - i) . = a 2 b 2 t [149] 

The solution of [1473/ meeting the initial conditions/ is 
6 = \ V 1 " cos [150] 

where 

% = a ^ s [151] 

is the nutational frequency. Comparison with [89] shows 
that the induced nutation has the same frequency as the free-body 
nutation. Since its amplitude is so small, it is expected to have 
a negligible effect compared with the free-body effect. 


The precessional angular velocity can now be obtained from [126] 


• cb . . crb6 

a = 7T- (u - u) - — — — 

_ -_2 o ■ sm y 
sin y 'c 


i 00 5, C1 - cos V’ 

sin Y 

1 A 


[152] 


The precession thus varies with time, ranging from zero to 
twice the average value . 


to 


2 ““s V sin \ 


3 Jr ,cr-l A 

2 a ("ET 5 cos T o 
s 


[153] 


As seen from the oi- frame, the motion thus consists of a steady 
coning of the spin axis around the orbit normal at the angular 


velocity (153] , with a very small superimposed circular motion, 

t 

with an amplitude ^'6 . 

In order to find the effects on the experiment of precessional 
motion, we neglect the nutation. Then [128] gives 

' * = “s [154] 

If the accelerometer sensitive axis is along the 2-axis, the 
test mass experiences a centrifugal acceleration given by 

a c = ~ 4 X ~ w p x sin2 Y 0 sin 2 u s t [155] 

from [121] . 


This is similar to the form [93] found for free-body nutation, 
except, that A is replaced by tocsin y Q and o) n by u s . One may then 
carry through a similar investigation of the stability of the 
solutions of the resulting Mathieu equation of motion. Analogously 
to [99] and o [98] , the coefficients in the standard form of the 
Mathieu equation are found to be in this case 


a = (1- - ) 
v 0) J 
s 



[156] 


and 



w sm y 0 
rP o, 2 

L 2r J 


[fc 2 ^) sin 2 y o ] 2 (J2/u 3 ) 4 


[157] 


By elimination of (ft/ai ) , one can obtain an equation showing 
the path followed by the operating point in the a-b plane as this 
parameter is varied. It is not difficult to see that this path 
is confined to the stable region near the origin in Fig. IV. 


It thus appears that gravity-gradient induced nutation and 
precession should not have a significant effect on operation of 
the experiment. The effects can, of course, be minimized by making 
Y 0 as small as possible, but stringent tolerances in the alignment 
of the spin axis with the orbit normal do not appear to be required 

A more detailed analysis nevertheless should be carried out, 
retaining higher frequency terms in the gravitational torque [112] , 
which were eliminated in the averaging process to obtain [119]. 

III. 3 C onclusions from Systems Analysis 

We have now considered in some detail the first four of the 
possible disturbances listed in Sec. III. 2. No effects have been 
found which might prevent operation of the experiment down to the 
design sensitivity of one part in 10 14 in the measurement of the 
Eotvos ratio. 

It appears to be possible to carry out a useful experiment 
in the payload bay of the shuttle, if it is feasible to locate the 
experiment reasonably close to the CM of the shuttle and if the 
RCS system can be used to fly the shuttle so as to follow the 
experiment. An average fuel flow of order 5 kg/hour would probably 
be sufficient for this task . • 

The STS (shuttle) provides an ideal vehicle for carrying out 
this experiment, which can make good use of the capabilities of 
the crewmen. If flown unmanned, the very sensitive accelerometers 
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would require caging to withstand the boost environment and some 
form of orbital gyrocompassing would be needed in order to erect 
the wheel spin axis to the orbit normal. Fui'thermore , if it is 
desired to measure the Eotvos ratio of several materials, a mechanism 

would be required for changing the test masses. Finally, experience 

. 32 

with the low-level accelerometers m orbit suggests that 

difficulties may arise which are very difficult to solve in an 

automatic system. The crewman may make final design choices as 

the result of on-orbit experience (e.g., in selecting servo compensation 

networks) , deploy the system in its optimum configuration, monitor 

the performance and investigate any anomalies observed, and modify 

the experiment (e.g., by installing new test masses) to extend the 

investigation. At the same time, the equipment may be considerably 

simpler than would be needed for an unmanned launch. 

© 

The development program for this experiment might involve three 
separate shuttle-born phases. In the first, component parts of the 
apparatus (particularly the accelerometer) would be tested in free 
fall, so that commitment to final development could be postponed 
until the design performance had been demonstrated. In the second 
phase, a limited experiment would be carried out in the payload bay, 
if further analysis shows that this is feasible and that significant 
simplifications could be effected by operation in this mode. Finally, 
(and particularly if any Eotvos anomalies were discovered in the 
second phase) , a fully autonomous experiment would be launched from 
the shuttle, at the maximum sensitivity which may .prove feasible. 
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CHAPTER FOUR 

DESIGN OF A MAGNETIC MICROBALANCE 
FOR ACCELEROMETER TEST 

IV. 1 Introduction 

As discussed in Chapter II, the heart of the orbital Eotvos 
experiment is a very sensitive dual accelerometer, the two test 
masses being constructed from the materials whose ratio of gravitational 
to inert mass is being compared. In low orbit and in the absence 
of external disturbances, the sensitivity of this accelerometer, 
in gees, is equal to the limiting accuracy in measurement of the 
Eotvos ratio. For example, to reach an accuracy of ^=10 in 
the experiment, the dual accelerometer must be capable of detecting 
a relative acceleration between the test masses of 10“^g. This 
is many orders of magnitude beyond the sensitivity of conventional 
accelerometers, so that design of the accelerometer system becomes 
the most critical task in development of the experiment. The 
frequency at which the acceleration must be measured can however 
be chosen by fixing the angular velocity of the apparatus, a 
freedom which can be used to discriminate against some sources of 
measurement noise . 


At the present time, both cryomagnetic and electrostatic 
techniques for suspension of the test masses are under study, the 
first primarily at Stanford and the second at M.I.T. 

A fundamental difficulty in the development of any sensitive 
accelerometer for use in space is that, if the device is to be 
tested on Earth, much stronger suspension forces than are needed 
in orbit must be provided. This problem has, in the past, generally 
led either to design compromises or to a reliance on analytical 
design techniques, without specific test of the final configuration 
before commitment to flight. 

In the Bell Miniature Electrostatic Accelerometer (MESA), for 
example, the gap between the test mass and the supporting electrode 

structure is comparable to that normally found in electrostatically 

© 

supported instruments for terrestrial application (i.e., thousandths 
of an inch) , so that it can be operated on Earth by applying not 
unreasonably high voltages. For operation in space, the electrode 
voltage is merely reduced, so as to provide support against the 
expected small transverse accelerations. The cross-coupling between 
forces along the suspension axes and the sensitive axis, which is 
due primarily to geometrical imperfections of the test mass and 
electrodes, electrode edge effects, etc., is thus no better than in 
terrestrial instruments, and little advantage is taken of the benign 
orbital environment to maximize performance. 


Most satellites designed for "pure gravity orbit" are examples 
of the other approach to this problem. Such satellites are driven 
so as to follow an unsupported test mass, which is protected from 
external forces . 22 This system can be regarded as an accelerometer 
which is sensitive to accelerations in all three axes. It is 
essential to minimize the forces (primarily electrostatic and 
gravitational) which are applied to the test mass by the satellite 
itself. This is usually achieved by careful analytical design, 
but there is generally little opportunity to test the resulting 
concept except by -an actual flight in space. 

An interesting approach to testing an electrostatically- 

33 

supported space instrument is that of de Bra et al , in which a 
test mass of sufficiently low density is used so that it may be 
floated in a heavy gas (SFg). This technique is not however 
considered applicable to the problem of testing accelerometer 
designs for the Eotvos experiment: it obviously will not work 

I 

at the low temperatures required in a cryomagnetic suspension and, 
at the higher temperatures which may be used with an electrostatic 
suspension, it introduces convective disturbances and damping forces 

In order to allow experimental investigation of design con- 
cepts for an electrostatically-supported accelerometer for the 
orbital Eotvos experiment, it would be highly desirable to simu- 
late free fall in the laboratory. The feasibility of using mag- 
netic forces for this purpose is investigated in this chapter. 
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Magnetic suspensions have of course been employed for a wide 
Variety of purposes. The present concept is unusual in that it 
is desired, not to provide stable suspension of the test mass, but 
to buck out its weight, using forces generated by appropriate sus- 
pension coils, in such a way that additional magnetic forces do 
not result when the test mass is displaced from its nominal posi^ 
tion. With the magnetic suspension providing the main support force 
and neutral stability, 'an electrode structure can be built to con- 
trol the test mass whose dimensions and voltages are similar to 
those which might be used in space . 

It is assumed here initially that the magnetic material em- 
bedded in the test mass is linear — - i.e., the magnetization vector 
induced in it is proportional to the external field. The reasons 
for this Condition and means for implementing it are discussed in 
Section IV. 5 below. 

■ - ■ o 

IV. 2 Basic Relationships 

The force on a sample of linear magnetic material in an ex- 

34 

ternal field B is given by 

F =-- k VB 2 1158] 

where k is a constant depending on the size, shape and 
magnetic properties of the object. If the body is diamagnetic 
(e.g., if it is superconducting), k is negative. 

The external field is given in terms of the vector potential A by 


i I I i l 1 . I 
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B = V x A [159] 

and the vector potential is given by an integral over 
the current density J in the coils: 

** ' k “ | d r ’ [160] 

where h° is the permeability of free space and r, r’ are 
vectors from the origin of coordinates to the observation point 
and to a source point in the coils, respectively. 

For constructional simplicity, it is proposed to build the 
suspension from a set of cylindrical coils of rectangular cross- 
section, with the symmetry axis (the z-axis) vertical. Figure VI 
shows a vertical section through one of the coils: the centroid 

of the cross-section is at radius b and height a above the origin, 
and the cross-sectional area measures 2(3 by 2a. The current density 
has only a circumferential component, which is 

j = m b “P < r ' < b+B 

■ a_ct < z' < a+3 - [161] 

~ 0 elsewhere 

where I is the current and n the number of turns/unit 

area. 

Because of the cylindrical symmetry, we may choose the observa- 
tion point in the r-z plane. From the geometry of the perspective 
sketch of the coil shown in Figure VII, we may write 
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• f < r ' z > = k ~ 1 ” 1 

= [(z-z 1 ) 2 + r 2 A- r' 2 - 2rr'cos £l€2] 

where <j> ' is the angle between the horizontal projections 
of r and r ' . 

To compute the integral [160] , we must write J in rectangular 
components , which are 


J x “ ’ J <f, sAn 


= J, cos <i' 

y <p r 


[163] 


It is clear from [1603 that A also has only a circumferential 
component, so that, at' the observation point, the integral over 
J x must vanish. The vector potential from all the coils can then 
be written ° 


0 b+0 a+a 27T 

A , = \ nl J r’dr’ J dz’ J f(r,z)cos t})’ dcf) 1 

a-a o 


4> 4tt r. ,. a 
Y coils b-p 


1164] 


The <p ' -integral may be expressed in elliptic integrals, 
but, for present purposes, it is more convenient to expand £(r,z) 
in a Taylor series in r about the symmetry axis, so that the 
integral becomes 


27T 


$ = J [f° ■+ rf (0 ,z) + — r f {0, z) + — x f (0,z) 
o 

. + ^r 4 f iv (0,z) + ■|- i r 5 f V (0,z) + — ] cos cp’ dtp' 


[165] 


■n 3 n F 

where f (0,z) = 


, so that 


r=o 


i 


( 


i 



f° = f(0,z) = [(z-z 1 ) 2 + r ' 2 ] ~V2 


[166] 


A straightforward computation of the indicated derivatives 
of [162] with respect to r yields 


f X (O f z) 

- f° 3 r'cos cj)' 


[167a] 

f 11 (0, z) 

5 2 2 

= 3f° r’ cos V' - f° 

3 

[167b] 

f 111 (0, z) 

7 3 3 

= 3[5f° r' cos (j)' 

3f° 3 r'cos <f> ' ] 

[167c] 

f 1V (0, z) 

9 .4 4,, 

= 3 [35f ° r cos <j>' - 

30f° 7 r' 2 cos 2 c|)' + 3f° 5 ] 

[167d] 

f V ( 0 , z ) 

= 15 [63f °' L1 r ■ 5 cos 5 (j) 1 

- 70f°^r 1 3 cos 3 <j> ’ + lSf^r'ccstj)'] 

[167e] 


When these expressions are used in [165] , all odd powers of 
cos integrate to zero. We are left with 


$ s TT[f° 3 r'r + f° 7 r' 3 - f° 5 r')r 3 

, 15.21 11 ,5 7 9 ,3 L _ 0 7 5 . 

+ — (— f° r' - — f r' + f° r')r + . . .] 
o o Z 

The vector potential may thus be written 


[168] 


\ = Q(z)r + K(z)r 3 + N(z)r 5 + ... 


[169] 


where 
_ P‘ 


su) - f t „i r r .^3 

4 L nj - j Jr' f° dr'dz' 
coils b-$ a-a 


[ 170 ] 


NL I 


b+8 a+a 


| W I nl / / [| f° 7 r' 4 - f° 5 r' 2 ]dr'dz- 


8 ^ *• » * l 4 

coils b-3 a-a 


[171] 


, c b+3 a+a 

/ lff U r. 6 

coils b~3 a-a 


q- f ° 9 r ' 4 + f° 7 r' 2 ]dr'dz' [172] 


If we differentiate [170] under the integral twice with 
respect to z (denoted ,by primes) , we obtain the important result 


Q" (z) = - 8 Kiz) 


[173] 


From [159], in cylindrical coordinates, the axial and radial 
components of the magnetic field are given by 

B - - -r— (rA , ) 

z r 9r 4> . 


2Q(z) + 4K(z)r 2 + 6N(z)r 4 + 


[174] 


B = - -K- A . 

r 9z 4> 


[Q'(z)r + K'(z)r 3 + N'(z)r 5 + ...] 


[175] 


Using [173] , the square of the magnetic field is then 


2 2 2 
B = B + B 
z r 


= 4Q 2 (z) + (Q' 2 (z) + 16K(z)Q(z))r 


+ 2(8K 2 (z) + Q'(z)K' (z) + 12Q(z)N(z))r 4 + 0(r 6 ) 


= 4Q 2 (Z) + (Q ' 2 ( z) - 2Q(z)Q"(Z))r 2 


+ 2(8K 2 (z) + Q' (z)K' (z) + 12Q(z)N(z))r 4 + 0 (r 6 ) 
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Since there is no term linear in r, there is, at the origin 
of coordinates, only an axial component of force (as one would 
expect from symmetry)', which is, from [158],* 


F = 8kQ°Q° ' 
z 


[177] 


where = Q(0), Q° ' = ~Q{z) 


dz x '“' z = o 


The objective of the suspension is to simulate free fall: 
the support force should be as nearly constant as possible (and 
equal to the weight of the suspended object) in a region around 
the origin. In Cartesian tensor notation, the force at a position 
Xj_ (i = 1,2,3, with the 3-axis vertical) is given by a Taylor 
expansion around the origin 


% = f? + s..x. + F . . x . 

1 1 1 D D iDk j 


■\ + 


[178] 


where 


F i = k !~' B2 i 

1 ox . x . =o 

1 1 


[179] 


s. . = k b 2 | 

3-D ox ;9x . X . ,X .=0 

1 D 3- 3 


[180] 


r - -v 
3-Dk 


2 3x. 3x.3x ‘x. ,x. ,x =0 

3- 3 k x j k 


[181] 


The simulation will be successful at the origin .if we can make 
the stiffnesses zero, and the "playing area" (i.e., the region 


♦Notice that, to give an upward force, Q° and Q 0 ' must be of opposite 
sign if k < 0 and of the same sign if k ) 0. 
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in which the simulation is acceptable) will be maximized if we can 
also make the second-order stiffnesses F . zero. In calculating 

1 JK 

these quantities, we take derivatives of [176] and then place r = 0, 
so it is clear that we can neglect the term in r^. In the present 
notation, [176] thus becomes 


B 2 = 4Q 2 + (Q ' 2 - 2QQ") (x 2 + x 2 ) + . 


[182] 


where Q is a function of x^ alone. 


The force is now 


F. 

i 


= k 


' 2 (Q‘ 2 - 2QQ" ) 

2(Q’ 2 - 2QQ")x 2 ^ 
8QQ’ - 2 QQ"' (x 2 + x 2 ) 


[183] 


which, at the origin, of. course reduces to [177] . 


The matrix of second derivatives of B 2 is 


9^B 2 

ax^Xj 


2(Q ,2 -2QQ”) 


2(Q* 2 -2QQ") 


-4QQ"' x 3 
-4QQ"* x. 


“4QQ'" X x -4QQ"' x 2 8 (Q' 2 +QQ") -2 (Q'Q'" +QQ ,m Cx^+x*) 


2 2 , 


[184] 


At the origin, the only surviving components are the radial 
stiffness 


S r S 11 " S 22 “ 2k(Q°’ 2 - 2Q°Q°") 


[185] 


i 
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and the vertical stiffness 


S z = S 33 = 8k (2°' + Q°Q°") 


[186] 


It is clear from these expressions that S z and S r cannot both 
be zero. Moreover, in conformity with Earnshaiw’s theorem, it is 
impossible to make the suspension stable in all directions if k is 
positive. We therefore choose to make the radial stiffness zero, 
relegating control along the vertical axis to an external stabiliza- 
tion system (see below) . The design goal is thus 

[187] 


q °» 2 = 2 Q°Q°" 


Note that, when this condition holds, the vertical stiffness 
becomes 

,2 


S - 12kQ° 
z 


[188] 


so that the system is vertically stable for diamagnetic test 
masses (k, <C 0) and unstable for k 0. A useful figure of merit 
for suspensions of this type is the characteristic length 


z = f /s 

z z 


1 0 ° 

3 Q° ' 


[189] 


which will be negative if k < 0 and positive if k ^>0. 

For a given support force, Z is maximized (i.e. , the vertical 
inhomogeneity is minimized) by making Q° as large as possible and 
Q° 1 small. Under the condition [187], 

e°" - f q°' 2 /q° = ] [190] 

and so must be small but finite. 


I 
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The second-order stiffnesses , calculated by taking the 
derivatives of [184] , form a three-dimensional matrix, but most of 
the terms vanish at the origin. By inspection, we find in fact 
that the finite terms are 

r il3 = r !31 = r 223 = r 232 = ^ll = F 322 = 

[1911 

r 333 = 4k(3Q o, .Q 0 " +‘ qoqohi > 

The expansion [178] of the force may now be written 
in vector form as 


F = k[8Q°Q° ' + 8 (Q° ' 2 +Q°Q° ") z - 2Q 0 Q 0 '" (x 2 +y 2 ) +4.(3Q° ' Q^'+SOQ 0 ’" ) z 2 ] i 


+ k[2(Q < *' 2 -2Q°Q 0, ')x - 4Q 0 Q°"' xy)i^ 


+ k[2 (Q° ' “2Q°Q°") y - 4Q 0 Q <>,,, xy)^ 


[192] 


where i ,i r i are unit vectors along the axes. 

—x — y — z 

If the condition [187] holds and if, in addition, we can 
achieve 


)Oflf __ 


[193] 


then the force reduces to 
F ~ k[8Q°Q 01 ' + 12Q° ' 2 z + 12Q 0 ’Q°"z 2 ]x 


~z 


[194] 


or 
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i - rn f + |(f) 2 i 


[195] 


Achievement of these conditions thus means that, to the 
second order in the displacement of the suspended object from the 
origin, no forces are generated in the horizontal plane. Further- 
more, the variation in the vertical force is independent of 
horizontal displacements, to this order. The quadratic term in 
[195] merely means that the vertical stabilization servo should 
have a slightly non-linear characteristic. 


The quantities needed for design of this system are obtained 
from [170] as 


. b+B a+a „ 

Q° = ~ y° l nl J f r' 2 ( z ' 2 +r' 2 )~ 3 / 2 dr'dz' 
coils b-3 a-a 


[196] 


■3 b+B a+a 

f I m / / r- 2 z -( Z ' 2 +r'V 5 / 2 dr-d Z - 

coils b-B a-a 


[197] 


n b+B a+a 

Q«" = fw l ml I r'hi^-r- 2 ) (,' 2 + r'V 7/; W'd,' 

coils b-B a-a 


[198] 


b+B a+a 


" ' ¥ l nI J J r' 2 z * (4z ,2 -3r' 2 ) (z • 2 +r • 2 ) ' 9/2 dr *d Z ’ [1 " ] 


coils b-B a-a 
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IV. 3 Ideal Coil Design 

To simplify selection of a configuration for the suspension, 
we assume at first that the coils are of infinitesimal cross-section. 
The above expressions then take on the simple forms 


Q° 

Q° 1 

Q°" 

Q °,., 




2 2 2 “ 3/2 

° I NI d (a +n ) 

coils * 

2 2 2 “ 5 / 2 
° l NI b a(a +b z ) 

coils 

2 22 22 " 7/2 

° S NI b (4a -b ) (a +b ) 

coils 

2 2222 ” 9/2 

° 2 NI b^'a (4a -3b z ) (a^+D ) 

coils 


[196a] 

[197a] 

[198a] 


[199a] 


where N is the number of turns in each coil. To reiterate, 

o 

the conditions we wish to achieve are [187] and [193] : 

Q° ' 2 = 2Q°Q°" [200] 


and 

Q° = 0 [201] 

At the same time, we wish to maximize the support force at the 
origin 

F° = 8kQ°Q° 1 [202] 

z 

and the length which characterizes the vertical homogeneity 
Z = |g°/Q° ' ;"[203] ' . 


IV. 3.1 Single-Coil Suspension 

Let us first consider the possibility of a suspension using 
a single coil. From the above expressions, the condition [200] 
gives, for this case 


a 2 - §b 2 


[204] 


There is thus, a point on the axis of a thin coil of radius 
f? 

b, at a distance ~ b from the plane of the coil, where the 

15 

transverse gradients of the magnetic force vanish. At this point, 
the characteristic length (203) is 


2 a 2 +b 2 _ 7 


a a 


b - .492 b 


[205] 


The support force is 


SkQ°Q° ’ 


_ 3 b 3 4 a 



,k]j ° 2 N 2 1 2 


3 4 ky o2 N 2 I 2 

. • ~ n 


•Tv T' 


0.247 


[ 206 ] 


It is clearly impossible to satisfy the condition [201] at 
this point. We therefore turn to suspensions containing two or 
more coils. 


V 


lV.3.2 Helmholtz Coils 

From [198a], Q°" vanishes is 


a = ib 
2 


[207] 


or if the coil lies on a cone with apex at the origin 


and semivertex angle 


4> h « tan"” 1 2 = 63.43' 


[208] 


Now suppose we have a pair of such coils, symmetrically 
located at Z = + a on the Helmholtz cone , with the same number of 
turns in each and the same current flowing. Then [196a] to [199a] 
give for the pair 


Q° = |u°NIb 2 (a 2 +b 2 ) 


2 y°Ni 
p?2 a 


[208a] 


QO » _ QO II _ QO III- Q 


[208b] 


This is, of course, the Helmholtz coil arrangement which is 
frequently used to produce a uniform magnetic field. 

If the current in the Helmholtz pair flows in opposite senses, 
we obtain 


= Q° " = 0 

? 7 77 “5/2 

= |*i°NI b^a (a z +b z ) 

„ v6 y°NI 

P72 


[210a] 


[210b] 


i 


I 


l 


J 


l 
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\ o « ii — 15 


y°NI b 2 a (4a 2 -3b 2 ) (a 2 +b 2 ) ” 9y/2 


= _ 48 y °NI 

5772 a 4 

Neither a Helmholtz nor anti-Helmholtz pair, by itself, can 
provide a lift force. 

IV . 3 . 3 Maxwell Coils 

From [199a], Q 0 "' vanishes if 


[210c] 


a 2 = rib ^ 
4 


[ 211 ] 


or if the coil lies on a cone with apex at the origin 
and semi-vertex angle 


(J) = tan 1 2 - 49.11° 

M /3 


[ 212 ] 


A pair of similar coils, symmetrically located at Z = + a o.. the 
Maxwell cone, with current in the same sense, yield 


Q° = iy°NI b 2 /(a 2 +b 2 ) 3//2 = 2 

2 7 3/2 a 


[213a] 


Q°" = 3 y°Nl b 2 (4a 2 -b 2 ) / (a 2 +b 2 ) 7/2 = 16 » 3 . 5/2 1L°NI [213b] 


7/2 


Ol — nOII! — 


Q = Q 


[213c] 


If the currents are in opposite senses, we obtain 


QO _ QO II _ QO Ml _ 


0 


[214a] 


[214b] 


Q° 1 = lu°NI b 2 a/(a 2 +b 2 ) 5/2 = 2(1) 5/2 


IV. 3. 4 Coil Combinations 

It should be clear that, by appropriate choice of Helmholtz 
and Maxwell coils, it is possible to provide independent control of 
Q°, Q° ' , Q° " , Q° , as follows: 

1) Q° alone is generated by a Helmholtz pair. 

2) Q°' alone is generated by an anti-Maxwell pair. 

3) Q°" alone may be generated by a Maxwell pair, together 
with a Helmholtz pair so connected as to subtract the 
Q° generated by the Maxwell pair. 

4) Q ° * n will be zero, as desired, if the suspension is made 
up of combinations of the above coil sets . 


Suppose we wish to generate Q°" alone. For a compact design, 

it will be convenient to use a Maxwell pair and a Helmholtz pair 

with the same radius b. For stability, we assume the same current 

I is flowing in both sets of coils (i.e., they are wired in series) 

From [213b] and [211] the Maxwell pair, with N turns in each coil, 

M 

produces 

Q°" = 334 y°NM I [215] 

7//2 


but it also produces an unwanted 



4 ^ j0 Nm i 
7372 b 


[ 216 ] 
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This can be compensated by a Helmholtz pair of N H turns each, 
with current flowing in the opposite sense, to give 


Q° = - 


4 y °Ntj 1 
72 b 


[ 217 ] 


from [209a] . We thus require 

7 ^/2 53 

v n h c = 1 - 651 ° M 


[ 218 ] 


I V . 3 . 5 A Maxwell-Coil Suspension 

The Maxwell coil has the advantage that Q 0 "' = .0 automatically. 
If we now consider a pair of coils, symmetrically located on the 
Maxwell cone, with the same current I flowing but with turns 

in the upper coil and N 2 in the lower, [196a] to [199a] give 


j Vl°I (N 1 +N 2 )b 2 /(a 2 tb 2 ) 


2 , . 7 N y°I 

U2 ( *l +N 2 )- b~ 


[219] 


o 999 5/2 

Q 0 ' = ~ y°I(N 1 -N 2 )b^a/(a^+b^) 

— jjj- ijq ) V° i 
7 5/2 (N 1 N 2 } 1T 


[ 220 ] 


Q°" = | y°K.N 1 +N 2 )b 2 (4a 2 -b 2 )/(a 2 +b 2 ) 7/2 
>= 7772 (N 1 +N 2 ) TT 


[ 221 ] 


[ 222 ] 


0 
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The condition [200] is now 


9(N 1 -N 2 ) 2 = 8(N 1 +N 2 ) 2 

The solutions of this equation are easily seen to be 

3 + /8 


V N 2 = 


34 


[223] 


[224] 


3-/8 

or the inverse of this. Which coil is to have more turns is 
readily determined by considering the force at the origin, which 
is * 

F° = 8kQ e Q° ' 


192 /3 ,„ 2 
-4 (N 1 


o °2 t 2 

. 7 2 . , v I 

* b J 


[225] 


The lower coil must have the greatest number of turns in the 
test mass is diagmagnetic (k < Q) . Test masses of equal but op- 
posite susceptibility could be supported by merely inverting the 
suspension. 

o 

With [224] , the support force becomes 


i O — 


0.138 | k j y° 2 N 2 I 2 /b 3 


[226] 


where N is now the number of turns in the larger coil. 
Comparison with [206] shows that, for coils of the same radius 
and number of turns and the same current, the single coil suspension 
produces a considerably stronger force than this Maxwell suspension. 
However, the characteristic length for the Maxwell suspension is, 
from [219] and [220] 


Z * 


14 


* 1 +N 2 


12 /I N l" N 2 


b * .714b 


1227] 


. i . • — 


I 


1 


1 


.1 


J 


89. 

which is considerably better than the single-coil case. In 
other words, the vertical inhomogeneity is less pronounced. 

As it stands, this suspension nominally meets the design re- 
quirements. However, in order to allow for construction tolerances, 
it is proposed to add an extra set of Helmholtz and Maxwell coils, 
outside the main coils, connected as discussed in the previous 
Section so as to produce Q°" only. Current flowing in these trim 
coils should have no effect on the lift force nor on the charac- 
teristic length, but can be used to trim the radial stiffness of 
the suspension to zero. 

The resulting coil configuration is shown in Fig. VIII. Some 
modifications to the turns ratios are required when the finite 
cross-sections of the coils are taken into account (i.e. , when 
[196] to [199] are used instead of [196a] to [199a] ) . 

IV. 3. 6 Vertical Servo 

Varying the current in the coils will not affect the condi- 
tions [200] and [201], but it will change the vertical force. 

This effect can be used to compensate for the gradients in the 
vertical force indicated by [195] , using the measured vertical dis- 
placement of the proof mass to control the current. This feedback 
control is essential for stability if the .proof mass is not dia- 
magnetic. 

To simplify current stabilization it is proposed to provide 
an additional pair of coils (called the servo coils), having turns 
n^ and in the ratio [224] and situated on the Maxwell cone. 
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For simplicity, it is assumed here that the servo coils are wound 


with the main coils, although they might have to be separated if 


the inductive coupling between the two sets of coils is excessive 


(especially if superconducting coils are used) . An accurately 


stabilized current I then flows in the main coils, while a current 


i(z) flows in the servo coils, in response to the measured vertical 


displacement z of the proof mass. The vertical equation of motion 


rs now 


mz = k' (NI + ni(z) ) 2 [1 + | + | (|) + . . .] - mg [228] 


where, from [226] , 


0.138 k 


[229] 


and N, n are the numbers of turns in the larger of the main 


and servo coils, respectively. Since it is assumed that the larger 


coil is the lower (upper) one, according as the proof mass is dia- 


magnetic (paramagnetic), k* is always a positive quantity. Recall, 


however, that Z is negative for diamagnetic proof masses. 


If the main coil current is adjusted so as to balance the 


weight of the proof mass at the origin, then from [226] 


mg = k* (NI) 


[230] 


and [228] becomes 


T /- — r /Bi ( 2 ) * * / ni(z) x 2 1/1 z j 1 /Z\ j_ \ 

2 /g - !2(- ljr — ) + (“Ni— ) 1 (1 + Z + 3 ( z> + •••> 


z 1 ,z. z 
+ Z * 1 <f> + • 


[231] 


We now write 


i(z) = 


f < C 1 I + C 2 <l> > 


[232] 
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i i 


in [231] . Retaining only up to quadratic terms gives 

z/g = (2C 1 + 1) | + (2C X + + 2C 2 + j) (|) 2 + . . . = C 

if 



[233] 


[234] 


If, then, we make 


1 ^ 2 Z + 12 + . . • ) 

this servo control will make the first- and second-order 
vertical force gradients zero at the origin. 


[235] 


Under these conditions, the suspension will support the magnetic 
proof mass in a neutrally stable region near the origin: no addi- 
tional forces are produced, up to the second order in displacements 
in any direction from the origin. - 


The objective of this suspension is to buck out the weight of 
a test mass so that it may be mounted in a weak electrostatic sus- 
pension of the type which might be used in the Eotvos experiment 
in free fall. In testing such an accelerometer, it is expected 
that the sensitive axis will be horizontal. The displacement de- 
tector required for the electrostatic support in this accelerometer 
may provide the signal for the vertical magnetic servo. 


It is clear that this system provides a high degree of isola- 
tion to the proof mass from vibrations of the support structure 








(e.g./ microseisms), since the suspension force is invariant up to 
the second order in displacements relative to the coils . The sys- 
tem is however sensitive to relative vibrations of the component 
parts. For example, suppose that the displacement detector for the 
vertical magnetic servo is displaced vertically by a distance e, 
relative to the coils, while the proof mass is at the origin (z=0) . 
The signal to the vertical servo will then be - e and the current 
in the servo coil will be, by [235] 


1 = 


i N i Ie 


2 n x Z 


[236] 


to the first order. According to [228] , the force produced 
on the test mass will be approximately 


F - k' (N-jT) * e/Z 
= mg e/Z 


[237] 


from £230] . The system must therefore be sufficiently rigid 
so that relative vibrations of its component parts, in the frequency 
range of interest, are very small compared to the characteristic 
length Z . 


IV. 4 Disturbing Magnetic Fields 

The vertical component of any external, spatially uniform, 
constant magnetic field merely adds a constant term, presumably 
small, to Q(z) in [174]. This will make a small change in the 
vertical force at the origin, according to [177], and it will al- 
so change slightly the radial stiffnesses, an effect which may 
be compensated by current in the trim coils so as to maintain 
the condition [187] . 
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The horizontal component of any external field is somewhat 
more troublesome. Suppose that there is an external field B g 
along the x^-axis, so that the magnetic field within the suspen- 


sion is 


B = 


sin 0 + 
r e 


B r cos 0 


B. 


[238] 


where jzS is an azimuth angle and B , B are given by [175] , 

i z 

[174] , respectively. The square of the field is then 


B 2 = B l + + B 2 + 2B e B r sin ^ 


2 2 

« B^ + B“ 
o e 


2B 


[Q* (x 3 )x 1 + K 1 (x 3 ) (x^ 4- y-2^ x i + •*•] [239] 


from [175]. B 0 is the field in the absence of the external 
field. 

The force on the proof mass when it is at the origin is given 
by [179] . The first-order effect of the external field is thus to 
produce an additional force in the x-^-direction given by 


= - 2kB e Q°' 


1 Be ma 


from [177] and [174] . 


[240] 


It appears from this equation that the field of the Earth, in 
particular, will produce an acceleration of the proof mass in the 
north-south direction whose magnitude, depending on the value 
achieved for B Q , may be of order one milligee . in using the sus- 
pension for accelerometer test, it would of course- be possible to 
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orient the sensitive axis east-west (or, more generally, perpen- 
dicular to the horizontal component of any external field) . If 
it is desired to simulate free fall to higher accuracy, it would 
be necessary to buck out the external field, for example by a large 
Helmholtz coil pair, outside the suspension with its axis horizontal 
The maximum field produced by this compensating coil need only be 
of order 0.2 gauss, and its direction and magnitude could be adjus- 
ted so as to produce a null field at the origin with the suspension 
switched off, as measured by a magnetometer. 

It is important to avoid time-varying external magnetic fields, 
especially in the frequency range of interest for testing accelerom- 
eters. In principle, it would be possible to measure the variable 
component of the external magnetic field with a 3-axis magnetometer 
situated outside the suspension (and far enough away so that it is 
not swamped by the main suspension field) . The signal produced 
could be used to vary the main (or servo) coil current, the trim 
coil current, and external Helmholtz compensation coil currents in 
such a way as to null the forces on the proof mass, but it is 
highly desirable to avoid this additional complexity. 

Prom [240] , the 60 Hz field due to power wiring will produce 
a jitter in the position of the proof mass, but the amplitude will 
be only a few hundred Angstroms even if the field amplitude reaches 
one gauss. This effect is therefore negligible. 

In the calculations presented here, it has been assumed that 
the suspension is constructed of ideal coils, of infinitesimal 
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cross-section. It is not difficult to show that the requisite 
suspension conditions [187] and [193] can be achieved using coils 
of finite cross-section. Realistic calculations require other 
design decisions about the suspension (in particular, whether it 
is to be constructed of superconducting coils or not) and there- 
fore are not reproduced here. 

Disturbing magnetic fields can be produced by two other 
mechanisms: (i) errors and distortions in winding the coils, so 

that the turns density is not uniform and the shape is not identi- 
cal to that assumed in the calculation; and (ii) distortions of 
the suspension field due to magnetic materials in its vicinity. 

To some extent these disturbances, being constant, may be corrected 
using the trim and compensation coils, but it is clearly desirable, 
from [240], to keep the fractional distortion of the field as small 
as possible. The limits to the simulation of free fall imposed 
by manufacturing tolerances have not yet been thoroughly studied. 

It is, in any case, important to build the support structure for 
the suspension from non-magnetic materials. 

IV. 5 Proof Mass Design 

In the work presented here it has been assumed that the mag- 
netic material of the proof mass is linear — i.e., that its 

. ■ 34 

magnetic moment may be written 

M = X V H [241] 

— m — 

where V is the volume of magnetic material and x m is its 
magnetic susceptibility. Under these circumstances, the change 
in energy in the suspension when the proof mass is introduced, 
with the coil currents held constant, is 
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u = M.I^ [242] 

where B = y H is the field without the proof mass. The 
o o— 

force exerted by the magnetic field is then 

1 ] X m ? 

F = VU = i V (M.B ) = ™ — V VB^ [243] 

o 

from [241] . This is of the form [158] . 

In order to achieve strong support forces, the obvious choice 
is to use a ferromagnetic proof mass. In this case, however, 

[242] is only a rough approximation, especially when hysteretic 
effects are included: the energy may depend on the path by which 

the proof mass is introduced to the suspension. Furthermore, the 
force on a ferromagnetic proof mass will not generally be of the 
form [243], even if [242] may be taken as sufficiently accurate. 

The use of magnetically soft and hard materials must be considered 
separately. ' 

IV . 5 . 1 Magnetically Soft Proof Mass 

If the magnetic field is sufficiently weak, a soft ferromagnetic 
object will exhibit a linear relationship similar to [241] (neglect- 
ing hysteresis) . A reasonable support force can then be provided by 
using a strong gradient in the magnetic field at the origin. Un- 
fortunately, a weak field and a strong gradient leads, according to 
[189], to a small value of the characteristic length Z of the sus- 
pension. 

In order to build a suspension of adequate performance, it is 
almost certainly necessary to use a field at the origin which is 
strong enough so that the proof mass is saturated. The magnitude 
of the magnetic moment will then be fixed, but its direction will 
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be along trio External field, or nearly so. We may thus write 


M = M c i 3 [244] 

where M is the saturated magnetic moment and i_ is a unit 
vector along the original field. Then [243] gives 

F = j M VB [245] 

so that the force is proportional to the gradient of the 
magnitude of the field’, rather than to the gradient of the square 
of the field. The analysis given in this chapter then fails. It 
may still be possible to build a suspension with the desired prop- 
erties using a soft magnetic proof mass, but the analysis will be 
considerably more complicated. 

. .. ■ * ■ 

XV . 5 . 2 Magnetically Hard Proof Mass (Permanent Magnet ) 

If the proof mass contains a permanent magnet whose magnetic 

moment is M and whose dimensions are much smaller than the 

— o 

characteristic length Z in the field, it will experience a torque 
given by 

T = M Q x B [246] 

In the absence of restraining forces, the magnetic moment 
will execute an oscillation about the field direction at an angular 
■frequency given by ■ 

u 2 - M o B/ms 2 [247] 

.where m is the mass and £ the radius of gyration about an 
axis perpendicular to the magnetic axis of the proof mass. As- 
suming for the moment that the period of this oscillation is very 
much less than the time constants characteristic of rectilinear 
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motion of the proof mass in the suspension, and that there is some 
damping of the angular oscillation (e.g., by eddy currents), the 
magnetic moment may be taken as being along the field direction in 
calculating the force on the proof mass. Equation [244] then ap- 
plies and, at the origin, the force is given by [245] as 

F = - mg = | M q VB [248] 

According to [174] , the field and field gradient at the origin 
are given by 

B z = 2Q° 

8 B / 9z = 2Q 0 ' [249] 

so that [247] and [248] give 


2 _ 2 g B .. 2g Q° 3gZ 
“ '~2 VB' 2 Q° ' 2 


[250] 


The angular frequency of oscillation is then that of a pendulum 
2 

of length £ /3Z and will generally be quite high if the magnet is 
small, as assumed in this section. 


An important difference between the hard and soft cases is 
that, if the magnetic moment is aligned with the field, a permanent 
magnet must physically rotate when it is displaced horizontally 
from the origin. It is easily seen from [174] and [175] that the 
magnetic field at a radial distance r from the origin makes an angle 

« “ ' Tz l251! 

A; proof mass containing such a magnet must therefore exhibit 
a rocking motion when it oscillates horizontally, an undesirable 
defect in the simulation of free fall. 


T ” ! 
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A possible solution to this difficulty arises from the ob- 
servation that, unlike the soft case, it is not necessary to 
provide a magnetic field at the origin when a permanent magnet is 
used. Apart from second-order effects unaccounted for in [246] , 
the torque vanishes at the origin if the field is zero there — - 
i.e., by [174] and [175], if Q° - 0. In the use of the suspension 
for accelerometer test, the magnetic moment may then be kept 
vertical by relatively small torques exerted by the electrostatic 
suspension system. We then have 

M = M i [252] 

where i is a unit vector along the vertical. Equation 
[242] then becomes 


by 


U = \ MB = M [Q (z ) - i Q"(z) (x 2 +y 2 ) + ... ] 


[253] 


2 “o~z “o 4 

from [174] and [175] . The force at the origin is then given 


E- 7 U lx,y,z -0 ’ M o iz " 

and the stiffnesses at the origin are 
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[254] 


[255] 


0 0 

The only non-vanishing terms in the 3-dimensional second-order 
stiffness matrix are 


r ll3 _ r 131 


f = O® ♦» » 

r 333 y 


r 223 r 232 r 322 "* r 311 


1 n o 1 1 i 

J Q 


[ 256 ] 
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It thus appears that a suspension using a permanent magnet 
can be constructed by choosing a suspension coil configuration 
such that 

Q° = Q°" - Q 0,M = 0 [257] 

with Q° ' finite. It was shown in Sec. IV. 3.4 that an anti- 
Maxwell pair (i.e., a pair of identical coils, symmetrically 
located on the Maxwell cone, with the same current flowing in 
opposite directions through them) meets these specifications. 

Notice that, in this permanent-magnet suspension, it appears 
to be possible to make the first- and second-order force gradients 
all zero. By Earnshaw's theorem, it would clearly be necessary 
to provide some form of feedback stabilization, but the lowest 
order terms in the destabilizing force appear to be cubic in the 
displacement from the origin. The performance of the system may 
thus be superior to that using a linear magnetic proof mass: in 

particular, it may be less sensitive to the type of relative vi- 
brations of the components discussed in Sec. IV. 3. 6. 

This interesting possibility deserves further investigation. 
For a permanent magnet, the principal defect of [242] is that it 
assumes that the proof mass is a perfect magnetic dipole. In a 
more complete analysis, multipole moments would be expected to 
interact with higher-order gradients of the magnetic field. As 
it turns out, however, the desired configuration is one in which 
at least the next two gradients of the field vanish, according to 
[257], so that the dipole approximation may prove quite accurate. 
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IV. 5.3 Diamagnetic Proof Mass 

We now return to the question of constructing a proof mass 
of magnetically linear material, for use in a suspension of the 
type described in previous sections. 

, . 34 

With the impractical exception of liquid oxygen, the dia- 
magnetic or paramagnetic susceptibilities of ordinary simple 
magnetic materials are too small to allow their use in a suspen- 
sion employing reasonable (kilogauss) fields and field gradients. 
As is well known, however, a superconducting body is perfectly 
diamagnetic - i.e., its permeability is zero. 


One solution to the problem, which may be practical if the 
suspension coils are superconducting (so as to maximize the 
achievable field and eliminate problems with heat dissipation) , 
is to embed a superconducting element inside the proof mass, 
which will otherwise be constructed of nonmagnetic materials. To 
eliminate torques on the proof mass, the superconducting element 
may be spherical. Its magnetic moment is then given in Gaussian 
units by^ 

M = *" <T- V B [258] 

— Off — — 

where V is its volume. The force is then given in cgs units 


by 


- 16ir 

so that, in [158] 


V VB' 


k = 




v 


[259] 

[260] 


The principal difficulties expected in the use of a super- 
conducting proof mass are those that have been encountered in the 
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development of superconducting gyroscopes. These are: 

(i) Elimination of trapped flux. In order that the suspended 
superconducting sphere not exhibit a permanent magnetic moment 
and hence anomalous forces, great care must be taken during the 
transition to the superconducting state. To avoid trapping the 
field of the Earth, it may be necessary to provide a mechanism 
whereby the proof mass may be kept outside the suspension, in 
the normal state, while the suspension coils are cooled and 
supercurrents established. The proof mass would then be inserted 
part way into the suspension, to a region where there is a 
definite field gradient in all directions, and allowed to cool 
until transition occurs. Only then would it be moved to the cen- 
ter of the suspension. 

(ii) Cooling of the suspended proof mass. If the proof mass is 
suspended in vacuum, it can be cooled only by radiation, a very 
inefficient process at cryogenic temperatures. The total power 
dissipation in the proof mass (for example, due to eddy currents 

in the non-superconducting parts) must be of the order of nanowatts. 
This condition may be relaxed by retaining a pressure of several 
torr of helium around the proof mass, so as to provide a convec- 
tive mechanism for heat dissipation. 

An alternative possibility, avoiding the problems of a super- 
conducting proof mass, is to use a pseudo-diamagnetic test body 

3 6 

of the type which has been investigated by Wilk . In this con- 
cept, the magnetic element in the proof mass would consist of 
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three concentric, orthogonal coils- of small cross-section, with 
a 3-axis (Hall effect) magnetometer at their center. The three 
coils carry currents delivered by servo amplifiers so as to drive 
the magnetic field at the center to zero. The element then behaves 
very nearly like a perfectly diamagnetic body. 

An active proof mass of this type of course requires a power 
source for the coils and internal electronics. However, in 
Wilk's analysis, the power requirements were found to be only of 
order 15 mW per gram of lift force, so that the possibility exists 
of building a proof mass with an internal power dissipation below 
one watt. If so, the power could be supplied by a photovoltaic 
cell forming one end of the proof mass, illuminated by a light beam 
with a power of order 5 watts. 

At the present time, neither the superconducting nor pseudo- 
diamagnetic proof masses have been adequately studied. The latter 
type, in particular, requires careful analysis in order to optimize 
the internal coil design for this application. 
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CHAPTER FIVE 

CONCLUSIONS AND RECOMMENDATIONS 

V. 1 The Magnetic Microbalance 

The analysis presented in Chapter IV has shown interesting 
possibilities for simulation of free fall in the testing of ac- 
celerometers. Of particular promise are magnetic suspensions 
using either (i) a linear magnetic element consisting of a Wilk 
pseudo-diamagnet ; or (ii) a permanent magnet with, its axis 
vertical, embedded in the proof mass. The suspension coil design 
to eliminate first- and second-order force gradients at the 
origin is quite different in these two cases and further work is 
required before a choice can be made between them. 

The original motivation for this study of magnetic suspen- 
sions was that it is impossible to test in the laboratory electro- 
static accelerometers of the sensitivity required for the orbital 
Ebtvos experiment, if it is necessary to buck out the weight of 
the proof mass with electrostatic forces. For the specific pur- 
poses of the Eotvos experiment, it may now be more cost-effective 


to test candidate accelerometer designs on an early flight of 
the shuttle (STS) , rather than to construct a suspension of the 

type considered here. It is estimated that a shuttle accelerom- 
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eter test facility could be constructed as a "minilab" weighing 

less than 200 kg. The direct launch cost, according to the NASA 

pricing formula and not including allowance for the services of 

a payload specialist, may be less than $200,000. 

To allow resolution of this issue, it is recommended that 
studies of the magnetic microbalance be continued to the point 
where a choice can be made between the above two types, the de- 
sirability of using superconducting coils in the suspension can 
be evaluated, and a realistic cost estimate for an operating 
system can be prepared. A parallel study should be undertaken 
of an accelerometer test facility for flight in the STS, so that 
shuttle interfaces may be identified, crew involvement specified, 
and preliminary estimates made of RDT&E expenses. Comparison of 
the cost and expected performance of the terrestrial and orbital 
accelerometer test options will provide the data needed for the 
orbital Eotvos experiment development plan. 

In making this decision, it should be noted that there are 
applications of the magnetic microbalance other than development 
of the Eotvos experiment. These include: 

(1) Use of the apparatus as a facility for testing electrostatic 
accelerometer designs for other space applications, in guidance 
and control of low-thrust vehicles, control of the attitude 
dynamics of satellites, establishment of "pure gravity orbit" 
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satellites, etc. Once the system has been built, its use for 
these purposes will be much cheaper and more convenient than tests 
carried out on STS missions. 

(2) Construction of a sensitive long-period seismometer/tiltmeter 
based on the microbalance. 

(3) Use of two microbalances to support the proof masses in an 

39 

electromagneticaily coupled gravitational-wave antenna for use 
on Earth. The isolation from microseismic and other disturbances 
provided by the system is particularly important for this applica- 
tion. Experience gained with such a terrestrial instrument, espec- 
ially if it proved possible to detect gravitational radiation from 
one of the expected cosmic sources (e.g., the Crab pulsar) wovild 
be of great val^.e in the design of a larger and more sensitive 
antenna for deployment in orbit by the STS. 

V . 2 The Orbital Ebtvbs Experiment 

Ac discussed in Section III. 3, the present state of develop- 
ment of the experiment allows confidence that it will prove pos- 
sible to compare the ratios of passive to inert mass of various 

14 

materials in orbit with a sensitivity of at least one part in 10 
This will allow a preliminary evaluation of the passive mass of 
energy stored in the weak interaction. Further analysis is how- 
ever required to ensure that none of the possible disturbances 
listed in Section III. 2 (or others) will preclude operation at 
the design sensitivity. 

The following immediate tasks have been identified for further 
development of the experiment: 
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(1) Continued analysis of disturbing forces. 

(2) Establishment of a firm candidate design for a dual electro- 
static accelerometer for this application. Of particular impor- 
tance is a realistic servo design to allow evaluation of the dis- 
crimination which is feasible against the identified disturbances. 


(3) Development of a plan for testing the accelerometer design, 
either using the magnetic microbalance or as an experiment on an 
early STS flight, as discussed above. 


(4) More detailed analysis of the feasibility of an Eotvos experi- 
ment to be carried out in the payload bay of the shuttle, including 
identification of shuttle performance limitations, shuttle and crew 
interfaces, etc. If the performance suggested by preliminary 
analysis is confirmed, finalization of system design for this ex- 
periment. 

(5) Firm estimates of the performance improvement which may be 
achieved by use of an autonomous, free-flying experiment which is 
erected in and launched from the shuttle. If the expected perform- 
ance warrants, detailed system design for this experiment. 

( (6 ) Development of a firm program plan, including cost and schedule 
estimates, for carrying out the various phases of the Eotvos ex- 
periment which are discussed above. 
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